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Abstract. Using the Flare Genesis Experiment (FGE), a balloon-borne observatory with 
an 80-cm solar telescope, we observed the active region NOAA 8844 on January 25, 
2000 for several hours. FGE was equipped with a vector polarimeter and a tunable 
Fabry-Perot narrow-band filter. It recorded time series of filtergrams, vector 
magnetograms, and Dopplergrams at the Cai 6122.2 A line, and Ha filtergrams with a 
cadence between 2.5 and 7.5 minutes. At the time of the observations, NOAA 8844 was 
located at approximately 5°N 30°W. The region was rapidly growing during the 
observations; new magnetic flux was constantly emerging in three supergranules near its 
center. We describe in detail how the FGE data were analyzed and report on the structure 
and behavior of peculiar moving dipolar features (MDFs) observed in the active region. 
In longitudinal magnetograms, the MDFs appeared to be small dipoles in the emerging 
fields. The east- west orientation of their polarities was opposite that of the sunspots. The 
dipoles were oriented parallel to their direction of motion, which was in most cases 
towards the sunspots. Previously, dipolar moving magnetic features have only been 
observed flowing out from sunspots. Vector magnetograms show that the magnetic field 
of each MDF negative part was less inclined to the local horizontal than the ones of the 
positive part. We identify the MDFs as undulations, or stitches, where the emerging flux 
ropes are still tied to the photosphere. We present a U-loop model that can account for 
their unusual structure and behavior, and it shows how emerging flux can shed its 
entrained mass. 
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1. Introduction 

Drawing on high-resolution vector magnetograms and Dopplergrams from 
the Flare Genesis Experiment (FGE) (Bemasconi et al., 2000; 2001), we 
report here on dipolar moving magnetic features that appear when flux is 
emerging and coalescing to form sunspots. Although it is common to see 
unipolar magnetic features flowing into growing sunspots (Strous et al., 
1996; Strous and Zwaan, 1999), we believe this is the first time that 
moving dipolar features have been reported in a sunspot formation region. 
Hereafter we will call these features moving dipolar features or MDFs, to 
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distinguish them from the moving magnetic features (MMFs) previously 
observed in sunspot decay. 

Sheeley (1969) first reported the MMFs that flow outward from mature 
sunspots, and Harvey and Harvey (1973) first described their magnetic 
properties. MMFs are about 1000 km in diameter and move about 1 km/s 
from a spot radially across its moat (Brickhouse and LaBonte, 1988) to 
collect in a ring of magnetic flux fragments surrounding the spot. In maps 
of the longitudinal magnetic fields, MMFs appear as either single poles or 
dipoles (Shine and Title, 2001). They appear to be an important 
mechanism of sunspot decay. No detailed explanation of the MMF 
phenomenon is established, but most observations are consistent with the 
‘sea serpent’ cartoon drawn by Harvey and Harvey (1973). In their article 
they suppose that discrete magnetic flux ropes are drawn away from the 
spot in partially submerged folds through mass motions originating in the 
Evershed flows of the penumbra. 

Nye, Thomas and Cram (1984) observed an MMF having a 400 m/s 
downflow, and we will see that downflow is a common feature in the 
sunspot-building moving dipolar features described in the present paper. 
The magnetic field intensity in MMFs has not been well established 
because of their small size, but the vector magnetic field observations 
reported here clarify why this is so. 

The FGE vector magnetograms reveal that the fields in MDFs are 
mostly horizontal and that they probably are areas where emerging flux 
ropes are stitched into the photosphere. That is, they appear to be U-loops 
(Spruit, Title, and van Ballegooijen, 1987), not the 1 2-loops usually 
associated with flux emergence (Parker, 1994). 

This paper is based on observations obtained during the Flare Genesis 
Experiment flight, which is described in Section 2 along with a summary 
of our data analysis techniques (Section 3). Readers most interested in the 
observations of MDFs in the emerging active region, NOAA 8844, should 
turn immediately to Section 4. We present our summary and conclusions 
in Section 5. The interesting relationships linking the emerging flux with 
the arch filament system seen in the chromosphere and with the X-ray and 
EUV structure of the overlying corona will be treated in another article 
(Schmieder et cil . , 2002). Georgoulis et a,l. (2002) present a statistical 
analysis of the over 300 Ellerman bombs that we recorded in NOAA 8844. 
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2. Instrumentation 

The Flare Genesis observations were obtained on January 10-27, 2000, 
during the FGE’s second Antarctic flight (Bemasconi et al., 1999; 2000). 
The principal element of the FGE is an 80-cm Cassegrain telescope with 
an F/ 1.5 ultra-low-expansion glass primary mirror and a crystalline silicon 
secondary mirror. The silicon efficiently conducts the heat, which is 
intense on the secondary mirror, into a thermal sink. Silicon also has low 
thermal distortion, so we believe the thermal distortion of the optics during 
the flight was negligible. 

A helium-filled balloon, provided by the NASA Scientific Ballooning 
Facility, carried the FGE to an altitude of 37 km where it could view the 
sun 24 hours/day without clouds or image distortion by atmospheric 
turbulence. Bemasconi et al. (2000) describe the focal plane 
instrumentation in detail, but in brief, it consisted of a polarization 
analyzer unit (two liquid-crystal variable retarders and a linear polarizing 
filter), a selection of 1.25 A order-isolation filters, a 0.16 A tunable Fabry- 
Perot etalon filter (Rust, 1994), refractive optics to present a 160" 0 
section of the Sun to a 1024 x 1024-pixel CCD electronic camera. The 
plate scale was 0.091" per pixel, and the image quality achieved was 0.5". 
We made observations of all four Stokes parameters in the red and blue 
wings of the Cal line at 6122.2 A. Unpolarized observations were 
obtained at 0.8 A in the blue wing of the Ha line at 6562.8 A. The Ca I 
observations yielded vector magnetograms, Dopplergrams and 
photospheric images. The Ha - 0.8 A images recorded activity in the 
lower chromosphere and in the arch filament system in the corona. 


3. Data acquisition and reduction 

We selected the Cal 6122.2 A line to derive the vector magnetograms 
because the 0.16 A bandpass of our Fabry-Perot etalon provides the best 
signal-to-noise performance with such a broad line. However, with an 
effective Lande factor of only 1.75, the line is only moderately sensitive to 
Zeeman splitting. 

The Stokes parameters were obtained from sequentially recorded sets of 
12 filtergrams at the polarization states I + Q, I -Q, I + U, I -U, I + V, 
and / - V, at ± 0.08 A in the wings of the Ca I line. The integration time 
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for each filtergram was 30 ms. Approximately 130 s were needed to 
acquire and store the full set of vector polarimetric measurements (12 
images). Following two such sets, a series of 12 chromospheric 
filtergrams was recorded. For these, the Fabry-Perot passband was in the 
blue wing of Ha, - 0.8 A off-band, and the integration time was 125 ms. 
Each sequence of two Ca I vector polarimetric sets and Ha observations 
required about 7.5 minutes. The sequence ran repeatedly for about 3.5 
hours (from 15:50 UT to 19:16 UT on January 25, 2000). A total of 56 
vector polarimetric sets and 324 chromospheric filtergrams were recorded. 

3.1. Basic data reduction 

The raw images were first corrected for camera dark-current offsets and 
then flat-fielded to account for pixel-to-pixel gain variations introduced by 
the camera electronics and optical system imperfections. To improve the 
signal-to-noise ratio, we resized the images to 512 x 512 pixels with a 
scale of 0.18" per pixel, corresponding to a field of view of 92.2" x 92.2" 
(66.3 x 66.3 Mm). We then carefully co-aligned the images to account for 
image motions caused by the residual telescope pointing errors and drifts. 

The telescope telemetry did not provide a high-precision position in 
latitude and longitude of the region observed. We determined its precise 
location by comparing a longitudinal magnetogram acquired by FGE at 
16:44 UT on January 25, 2000 with a Kitt Peak synoptic magnetogram 
recorded between 16:21 UT and 17:16 UT of the same day. 

3.2. Vector magnetograms 

We obtained Stokes images, Q/I, U/I, and V/I, by calculating the 
normalized difference between measurements of opposite polarization 
states. For example: 

V/I=[(I+V)-{I-V)]/[(I+V) + {I-V)l (1) 

The magnetic field components longitudinal (B,,) and transversal (B ± ) 

to the line-of-sight direction (LOS) were calculated following Ronan, 
Mickey, and Orrall (1987): 

B„ = L (V//), B ± =Ttj(Q/l) 2 +(u/l) 2 , (2) 

where L and T are two constants that depend on the magnetic sensitivity of 
the spectral line used, on the exact location of the Fabry-Perot passband, 
and on the instrument polarimetric efficiency. We derived L and T by 
comparing the FGE data with vector magnetograms obtained by the 
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Imaging Vector Magnetograph (IVM) of Hawaii (Mickey et a!., 1996), 
which was observing the active region at the same time. The IVM has a 
lower spatial resolution than FGE but it records full Stokes profiles, and 
the longitudinal and transversal field components are derived by means of 
an inversion code that fits the Stokes profiles. The advantage of 
comparing the FGE vector polarimetric measurements with the IVM 
inverted vector magnetic data is that we do not risk overestimating or 
underestimating the effective magnetic field strength in either B„ or B ± . 

Before comparing the two data sets, we resized, rotated, and co-aligned 
the FGE vector magnetograms to match the IVM resolution and field of 
view. From the comparison we derived the following calibration factors: 
L — 26100 ± 450 Gauss, T = 4500 ± 100 Gauss. 

The magnetic field zenith angle % i.e. the vector field inclination with 
respect to the LOS, and azimuth % were calculated, again following 
Ronan, Mickey, and Orrall (1987), as follows: 

r = arctan(^J, ^larctan (3) 

Because of the symmetry properties of the transverse Zeeman effect, 
another possible solution for the azimuth angle is: Zi = Z+ 180°. This 
180° ambiguity is not easily resolved. Many techniques have been 
proposed (see Metcalf, 1994 and references therein, as well as Skumanich 
and Semel, 1996). However, most techniques work well only for active 
regions that do not deviate too much from a potential field configuration. 
Active region NOAA 8844 contained several areas where the field 
structure was highly complex and non-potential. A totally automated 
algorithm would not be able to properly resolve the 180° ambiguity for 
those areas. 

We adopted an approach that combines a potential field solution with 
an empirical estimate of the correct field orientation based on simple 
geometrical considerations. Out of the 56 recorded vector magnetograms, 
we selected one, which was kept as reference, and we resolved the 180° 
ambiguity for it. Initially, we compared the field vectors with a simple 
potential-field model derived from the longitudinal magnetogram. In this 
model we just considered the two main polarities in the active region and 
neglected the magnetic flux between them. We chose the field azimuths 
that minimized the angular difference with the model azimuths. The 
resulting vector magnetogram was subsequently filtered with a 
“smoothing algorithm” that checks for maximum continuity between 
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adjacent vectors to create a smooth azimuthal distribution of vectors 
across the entire region. We then visually inspected the azimuth map and 
corrected manually any local geometrical inconsistencies. The vector field 
map was then transformed from the LOS reference system to the local, 
heliographic reference system. In the local reference system the zenith 
angle is expressed with respect to the local vertical and the azimuth with 
respect to the local N-S/E-W directions. The transformation into the 
heliographic reference system was helpful in revealing vectors that were 
incorrectly chosen due to the = 30° off vertical viewing angle of the region 
in the LOS system. See Gary and Hagyard (1990) for a discussion on the 
effects of perspective. The transformed vector magnetogram was again 
visually inspected and corrected for local geometrical inconsistencies. We 
corrected the rest of the vector magnetograms sequentially by comparison 
with the preceding vector magnetograms, starting with the reference one. 
We assumed that the vector directions changed less than 90° during the 
times between consecutive magnetograms. After the sequential 
comparison procedure, a few vector magnetograms still needed some 
manual correction to account for isolated local errors. 

Despite all our efforts to solve the 180° ambiguity problem, we still can 
not be completely sure that we chose the right solution across the entire 
field of view and for all 56 frames. However we believe that at worst, the 
ambiguity was incorrectly resolved only in a few isolated areas. 

3.3. PHOTOSPHERIC DOPPLER MAPS 

We derived photospheric Dopplergrams by calculating the normalized 
difference between filtergrams recorded in the two wings of the Cal line. 
We averaged all the images recorded in each wing during each 12-image 
set and over two consecutive sets: 

D b = (< Qb +Qb+ U b + U b + V b + V b let 1 + (• • -)set 2 » (4) 

i>,=(3; + a-+c/;+t/; + F; + r,^ 1 + (...) K12 , (5) 

where the indices ‘b’ and ‘r’ mean blue and red wing, respectively. This 
was done to improve the signal-to-noise ratio, to remove some of the 5- 
min oscillations effects, and to cancel out image-to-image intensity 
variations introduced by the vector polarimetric measurements. The 
Doppler image was obtained as follows: 

D = (D b - D r ) / (D h + D t ). (6) 

A positive value in D corresponds to a red shift in the line. To further 
reduce/remove 5-min oscillation effects and to reduce image-to-image 
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noise, we applied a space-time Fourier filter. The filter combined a high- 
pass filter, for the noise, with a bandpass filter, for the intensity 
oscillations. The final time series consists of 28 Doppler images with a 
time step of approximately 7.5 minutes. 

Each Doppler image D is just a map of the intensity difference between 
the blue and the red wings of the Ca I line. To estimate the longitudinal 
velocities V in the photosphere, we assumed a linear relationship between 
D and Doppler velocity V, i.e.: V — c D, where c is a calibration constant. 
We derived c by comparison with the velocity maps provided by the 
inversion of the IVM data: c = 38 ± 3 km/s. 

3.4. Error analysis 

We performed a statistical analysis of the errors in the Stokes parameter 
measurements and derived the errors in the inferred magnetic fields. From 
analysis of the polarization images, we estimate the following rms noise in 
the polarization measurements: Stokes V/I : ±0.0023, Stokes Q/I : 
± 0.0018; Stokes U/I : ± 0.0018. For Stokes Q and TJ the rms noise values 
are slightly smaller because we applied a weighted running average to the 
Stokes images to increase the signal-to-noise ratio even though this 
averaging slightly reduced the spatial resolution for the linear polarization 
maps. The Stokes Q and U maps show signals almost a factor of two 
smaller than Stokes V, therefore the effective signal-to-noise ratio is still 
considerably larger than for the V maps. 

From the errors in the Stokes measurements, we estimate the systematic 
calibration error in longitudinal field at 70 G with 60 G rms noise and a 
threshold of 55 G. For the transverse field the estimated systematic error 
is 85 G with 60 G rms noise and a threshold of 180 G. Here, the rms noise 
is defined as the error when measuring the same quantity multiple times. 
The threshold is the smallest value that can be measured above the 
background zero level noise. 

The error in measuring the vector field direction is highly dependent 
upon the strength of the magnetic field. The stronger the field, the more 
accurate the direction estimate. From our calculations, it appears there is a 
maximum error at each field strength. But the average error is always 
between 0° and the maximum value. Table I shows the change in average 
and maximum errors for the zenith and azimuth angles as a function of 
field strength. In the table, \B\ is the absolute value of the vector field and 
7?trans is the vector field component transverse to the LOS. 
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TABLE I 

Error estimates in the direction of the magnetic field vectors 


Zenith Angle Errors 
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From the statistical distribution of the field strengths, we estimate that 
the typical \B\ observed was about 1000 G. Therefore, we can say that the 
typical average error in estimating the zenith angle is about 3°. The 
maximum error is about 6°. For the transverse components, the typical 
value was probably 400 G. Therefore, for the azimuth angle we have 
typically a 5° error in average and 8° maximum error. 

These are just the statistical errors derived only from known 
measurement errors. It is extremely difficult to evaluate the errors in the 
vector direction that are related to the specific type of observation we 
made, i.e. measuring the polarization at a single wavelength position with 
a single and rather broad bandpass filter. We do not consider magneto- 
optical effects, but we do not expect that they are important as far out in 
the wing of the line as we were. In addition, we cannot evaluate the 
effects on the magnetic measurements due to Doppler shifts of the spectral 
line and other effects that may change the line profile. A comparison of 
the magnetic vector fields derived from FGE data and the corresponding 
IVM vector fields reveal that the zenith angle differences are mostly 
smaller than 5-10°. The differences in the azimuth angle estimates are 
usually bigger but they are also generally below 10°, with some areas 
(especially in penumbral regions) where the difference can rise up to about 
20 °. 

We estimate the errors in the Doppler velocities of about ± 0.2 km/s 
with an rms noise of 0. 1 1 km/s. We want to emphasize that the derived 
Doppler velocity maps probably underestimate the true line-of-sight 
velocities. Only measurements at multiple wavelengths across the spectral 
line can provide highly accurate velocity values. 



MOVING DIPOLAR FEATURES IN AN EMERGING FLUX REGION 9 

4. Observations 
4.1. Flux emergence overview 

Active region NOAA 8844 had only two, small, stable sunspots until 
January 24. Rapid growth started late on the 24 th or early on the 25 th . The 
region continued growing for only about another 24 h after January 25. 
Our observations (Figure 1) began on the 25 th at 15:50 UT and ended at 
19:16 UT that day. At the time of the observations the active region was 
located at approximately 5°N 30° W. 

Figure 1 shows NOAA 8844 at 17:57 UT as it appeared in the wing of 
the Ca I line (l.A), in a longitudinal magnetogram (l.B), in a Dopplergram 
(1C), and in the blue wing of Ha (l.D). Note the thin dark fibrils in the 
center of the region, evident in the Ca I filtergram. These fibrils are 
characteristic of flux emergence in growing sunspot regions (Bray and 
Loughhead, 1964). The Dopplergrams show that the dark fibrils (each ~ 
3000 km long) were rising through the photosphere with typical velocities 
of at least 1 km/s. Flux emergence was episodic, with each dark fibril and 
its signature of upwelling in the Dopplergrams being visible for 5 to 10 
minutes. The Ha images show signatures of an arch filament system, a 
characteristic feature of growing active regions (Bruzek, 1967). The tops 
of the arches are prominent in the Ha sequence, because they were rising, 
and while we cannot measure the exact rate of growth, the fact that they 
are clearly seen at 0.8 A in the blue wing indicates an upward velocity of 
about 35 km/s. Also prominent in the Ha images are the Ellerman bombs 
(Ellerman, 1917; Kurokawa et al., 1982; Nindos and Zirin, 1998), which 
have very broad, bright Ha wings and in our filtergram appear as 
brightenings roughly 1.5 arcsec in diameter. Georgoulis et al. (2002) 
analyze the Ellerman bombs in NOAA 8844 in detail. 

There were three principal sites of flux emergence, as indicated by the 
large dashed ellipses in Figure 3. One (labeled with a ‘C* in Figure 3.1) 
was in a supergranule adjacent to the follower spot and the two others 
were in supergranules just east of the center of the active region. 
Although we use the word, ‘supergranules,’ we note that they were only 
about 10 000 km in diameter instead of the more usual 30 000 km. 
However, the velocities - upwelling at the center and downflow along the 
edges - are characteristic of supergranules. The horizontal velocity 
pattern, with 0.3 - 0.5 km/s flows from the center, is also typical of 
supergranules. The vertical velocity pattern was determined from the 
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Fig. 1. Four views of NOAA 8844 recorded by FGE at approximately 
17:57 UT on January 25, 2000. (A) Filtergram in the blue wing of the Ca I 
6122 A line. (B) Map of the line-of-sight longitudinal magnetogram. 
(C) Dopplergram; bright corresponds to a blue shift, i.e. a flow motion 
towards the observer. (D) Filtergram at 6562 A: Ha - 0.8 A. The images 
have a field of view of 90x90 arcsec and the tick mark separation is 2". 
The orientation of the AR and the Sun center direction are shown in the 
lower left comer of image A. 

Dopplergrams. Cross-correlations of the pattern of granulation visible in 
successive Cal line wing images determined the horizontal velocity 
pattern (Shine, 2001). A check of magnetic element motions vs. the 
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Fig. 2. (A) Color-coded map showing the orientation (azimuth) of the 
transverse magnetic field component. The color wheel in the lower left 
comer explains the vector direction (imagine the vector pointing from the 
center of the wheel towards the rim. The black areas indicate the location 
of the spot umbrae and the contours outline the penumbrae. The white 
areas do not show any data. (B) Same as A but with the field vector 
orientation expressed by arrows plotted over the corresponding 
longitudinal magnetogram. Tick mark separation is 2". 

granulation horizontal velocities showed general agreement, so the surface 
motions appear to have been largely determined by the convective flows 
of the supergranulation. 

To a first approximation, we can say that the magnetic fields emerged 
in a wide band of horizontal flux ropes between the two major spots. 
Figure 2, which shows the vector field orientations, reveals that the 
direction of the transverse field in the flux emergence region was 
remarkably uniform, as shown by the purple and red area in the figure. 
Purple corresponds to fields tilted ~ 20° counterclockwise from the line 
connecting the leader and follower spots and red corresponds to » 60° tilt. 
The strength of the transverse component in the flux emergence region did 
not vary greatly. 

In general, negative poles flowed toward, and merged with, the 
follower spot, and positive poles flowed toward, and merged with, the 
leader spot. Typical velocity of the flow was 0.3 km/s. This emerging 
flux region was similar to the one studied recently by Strous et al. (1996) 
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and Strous and Zwaan (1999) and to the three regions studied by Lites, 
Skumanich, and Martinez Pillet (1998). However, the high resolution of 
our observations and the unvarying image quality from frame to frame 
revealed an important new fact: many of the emergent magnetic features, 
as seen in the longitudinal magnetograms, appeared to be dipolar. 

4.2. Magnetic properties of moving dipolar features 

Over an interval of 206 min, we observed approximately 20 MDFs in an 
area of about 5 x 10 8 km 2 . They were pairs of magnetic poles, each with 
an area of about 5 x 10 5 km 2 , that moved together as obviously joined 
features. Most of the MDFs first appeared near the center of one of the 
three supergranules located west of the negative polarity, trailing sunspot 
(see Figure 3.1). Most migrated at a nearly constant velocity between 0.3 
and 0.8 km/s from the supergranule centers toward the edges and/or into 
the sunspot. However, for short periods an MDF could dash at almost 1 
km/s. We interpret the sudden accelerations and decelerations as 
interactions with the surrounding flow fields. This idea is supported by 
the observation that frequently, during the changes in speed, the MDFs 
also changed shape, mainly stretching or shr inkin g 

In the supergranule cell nearest the follower spot (ellipse labeled ‘C’ in 
Figure 3.1), only a few MDFs were moving westward, i.e. toward the edge 
of the cell opposite to the spot; the majority moved towards the spot and 
eventually merged with it. In all cases the polarity of the part closest to 
the spot was opposite to the polarity of the spot itself. This is opposite to 
the orientation in most of the dipolar MMFs observed in decaying 
sunspots (Harvey and Harvey, 1973). 

Figure 3 shows a time series of longitudinal magnetograms that follows 
the position and evolution of two MDFs (ellipses labeled ‘a’ and ‘b’) as 
they approached the negative-polarity trailing spot. In the figure it is 
possible to see several more MDFs than the two we highlighted. 

MDF ‘b’ approached the spot and merged with it. Once the MDF drew 
close to the spot, the positive feature nearest the spot elongated slightly 
then gradually disappeared (Figures 3.2 and 3.3). Conversely, when the 
negative part approached the spot, part of the spot itself appeared to 
extend towards it, to merge with it, and then to retract again (Figures 3.4 
and 3.5). When an MDF approached the west edge of the supergranule 
(pointed out by the white arrow in Figure 3.2) from the east, its negative 
part disappeared and its positive part remained. Similarly, MDFs 
approaching from the opposite side left only a trace of their negative 
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Fig. 3. Sequence of longitudinal magnetograms showing MDFs near the 
following (negative) spot of NOAA 8844 (black area at top). Most MDFs 
appeared in the supergranules outlined by the large dashed ellipses (frame 
1). The small ellipses, marked by ‘a’ and ‘b’ in frame 1, follow the 
evolution of two such MDFs while they migrate towards the spot. The 
dashed ellipse labeled with C in frame 1 is the supergranule cell discussed 
in the text. The arrow in frame 2 indicates the location of the western 
boundary of the supergranule C. 

fields. The movie 1 of the magnetic evolution at the supergranule cell 
boundary gives the impression that positive and negative fields were 
canceling, which could mean they were reconnecting and/or submerging. 

The lifetimes of MDFs varied from tens of minutes to hours. For 
example MDF ‘a’ in Figure 3 is the most persistent one in our recordings. 
It survived for almost the entire 3.5 h of the sequence. In the last frame 
(Figure 3.6) its positive polarity appeared to extend towards the negative 
spot and it was probably about to be swallowed by the sunspot. Most 
MDFs disappeared only when they reached the edge of a supergranule or a 
sunspot, but a few faded and disappeared before reaching such a 
boundary. Figure 4 shows five typical MDFs. Three of them (labeled 1, 


See web page: sd-www.jhuapl.edu/FlareGenesis/Science/results.html#Movies 
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Fig. 4. Longitudinal magnetograms with ellipses indicating the location of 
the five MDFs discussed in the text. The arrows show the direction of 
motion of the MDFs. The tick mark separation is 1 ". 


4, and 5) moved toward the follower spot, located on the upper edge of the 
images, and two moved away from it (2 and 3). 

Figure 5 shows the vector fields in the five MDFs. In each MDF, the 
vectors point to the left in the positive and negative parts, revealing that 
each MDF vector magnetic field is a shallow U-loop. Also in each MDF 
the vectors in the negative part are steeper than the vectors in the positive 
part. 

For each MDF we derived the dimensions and field characteristics 
along a line that goes through the center of its two magnetic parts. Along 
this line we calculated the average for four parameters over a width of 
seven pixels 1.2" « 1000 km). The four parameters are the absolute 
field strength \B\, the field strength component parallel to the local vertical 
Bn local, the local zenith angle yiocai, and the line-of-sight velocity. Their 
variation along the MDF length is shown in Figure 6. Each line type 
represents a different MDF. We aligned the traces so that the neutral lines 
at the centers of the MDFs coincide. The plot shows that all the MDFs 
had similar characteristics. Each MDF attained a maximum of the 
longitudinal field strength at around 200 G, and each reached the peak 
approximately 1" from the neutral line that bisected the respective MDF 
(Figure 6.B). In contrast, Tlocai showed a clear asymmetry with respect to 
the neutral line (Figure 6.C). The field lines in the negative part of each 
MDF were more perpendicular to the local surface than the ones in the 
positive part. The field angle peaked on average at 50° in the negative 
part, while its peak in the positive part was only 20°. 
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Fig. 5. View in perspective of the magnetic vector field geometry for the 
5 MDFs of Figure 4. In each case, the vectors (yellow arrows) point 
downward and to the left in the negative part and upward and to the left in 
the positive part. 
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Fig. 6. Properties of the five MDFs highlighted in Figures 4 and 5. (A) 
Absolute value of the vector magnetic field strength. (B) Field strength of 
the vector component parallel to the local vertical (longitudinal field 
strength in the local coordinate frame). (C) Field inclination (zenith angle) 
in the local reference system: -90° is vertical down, 0° is horizontal, 90° is 
vertical up. (D) Plasma velocity parallel to the line-of-sight direction. 
Positive is towards the observer (upflow). 
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In general, the absolute field strength appeared to be greater in the 
positive parts than in the negative parts (see Figure 6.A). Moreover, the 
increase in \B\ included even the fields just outside of each positive part. 
We believe that this is mainly due to the proximity of the positive parts of 
the MDFs to the sunspot and its stronger fields. This effect was more 
prominent in MDFs 4 and 5, since they were closest to the spot (compare 
their positions in Figure 4). On the other hand, MDF number 2, which 
was located very far from the spot, showed a quite flat field strength of 
about 220 G over its entire length. By just looking at |B| maps it is very 
difficult to identify the MDFs. This means that generally the magnetic 
field strength within the MDFs was not remarkably different from the field 
strength in the adjacent areas. What made the MDFs visible in the 
longitudinal field maps was their vector field inclination. In the 

surrounding areas, the field lines were nearly horizontal, on average —5° 
(i.e. slightly negative polarity), while the field line inclination within the 
MDFs reached relatively high angles with respect to the surface. We 
therefore believe that the MDFs are not single isolated dipoles emerging 
from below the surface, but most likely dips in the canopy of horizontal 
field lines that filled the entire area of emerging flux. 

4.3. Material flow in MDFs 

An important feature of the MDFs is their Doppler signature, which is 
shown in Figure 6.D. All five of the MDFs there had a net downflow with 
respect to the surroundings. In particular, the Doppler curves show two 
clear dips coincident with the location of maximum longitudinal field 
strength and also the highest field angle with respect to the surface. This 
could be evidence for drainage along the field lines of the dense plasma 
that was carried into the upper photosphere / lower chromosphere during 
flux emergence. Note that the downflow was not confined to the parts of 
the MDFs with the strongest longitudinal field. Each whole MDF was 
sinking, according to the Dopplergrams. 

We have already remarked on how the MDFs interacted with the 
negative sunspot, namely, the positive parts disappeared and the negative 
parts merged with the negative fields of the spot penumbra. After each 
merger, we could not distinguish the MDFs in the magnetograms, but in 
the Dopplergrams the downflow signature persisted, and one can follow 
the MDFs as they pushed into the penumbra even though the penumbra 
had a strong outward Evershed flow. We are reminded of Parker’s article 
(Parker, 1978) predicting that draining the material of magnetic flux ropes 
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Fig. 7. Enlarged images of the Dopplergram (A), contrast enhanced Ca I 
line wing filtergram (B), and LOS longitudinal magnetogram (C) of the 
western edge of the trailing sunspot recorded at 16:07 UT. The ellipses 
indicate the location of two MDFs moving towards the spot. They are 
linked by a dark, rising fibril whose location is pointed out by the arrow. 
The fibril is a dark thread in the filtergram (B) and is associated with an 
upflow visible in the Dopplergram (A). Figure 8 sketches the field lines 
that could produce this particular configuration. 

will initiate enhancement of sunspot field strength by convective collapse. 
We have not yet examined our data for magnetic field enhancements 
precisely associated with MDF arrival at the penumbra’s edge, but we did 
note that the photosphere near the spot darkened as the MDFs arrived. 
Our results are consistent with Lites, Skumanich, and Martinez Pillet 
(1998) who found downflows on edges of sunspots in an emerging flux 
region. 

Figure 7 illustrates the relationships linking the dark threads of 
emerging flux, the upward flows, and the MDFs. In the figure we can see 
two MDFs that appear to be connected by a dark thread in the 
photospheric image (indicated by the arrow in Figure 7.B), co-located with 
a white thread in the Dopplergram (Figure 7. A) indicating an upflow. The 
MDF near the center of Figure 7.C is about to merge with the spot, and its 
positive part is stretching. Note in Figure 7.A the two dark patches 
(downflows) that are exactly co-located with the two more vertical 
magnetic parts of the lower MDF. The figure shows what we found to be 
the general case: MDFs with their entrained downflows appeared at the far 
ends of the dark threads in the continuum pictures. 
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Fig. 8. Sketch summarizing the relationship between emerging dark fibrils 
and the MDFs seen in longitudinal magnetograms. The curve B indicates 
the magnetic flux rope. 


5. Summary and Conclusions 

A summary of our findings: moving dipolar features have been discovered 
in emerging magnetic flux. They move horizontally toward sunspots and 
supergranule boundaries. All the MDF pairs, wherever they occurred in 
the emerging flux region, were oriented opposite to the main active region 
polarities, i.e., each positive (leading polarity) part was eastward, not 
westward, of its corresponding negative (following polarity) part. The 
MDFs were remarkably similar in size, field strength, and orientation. 
Downflows were persistent at each pole. The MDFs were bom and are 
sustained in regions with nearly horizontal fields (average elevation « 5°). 
They expired, apparently with a downward plunge, at the edges of 
sunspots and supergranules. Vector magnetograms (Figure 5) show that 
the MDFs are really magnetic U-loops stitched into the upper photosphere, 
where the Ca I line is formed. 

Figure 8 summarizes our interpretation of the data. We emphasize that 
the principal component of magnetic field in MDFs is the horizontal field. 
The observations show that bulges in nearly horizontal flux ropes at the 
upper photosphere create dark threads in the white-light images. The 
rising segments of flux rope are about 3000 km long. At a distance of 
another 3000 km, we find connected MDFs, whose poles are 
approximately 1000 km apart. Plasma in the dark threads would rise, cool 
(as indicated by the darkening), become more dense, and then slide 
downward into the MDFs. The fact that all the MDFs were similar in size 
and field orientation suggests that they result from a hydrodynamic 
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instability with a characteristic wavelength of approximately 1000 - 3000 
km. 

Perhaps the Kelvin-Helmholtz instability (KHI) (Min, 1997; Ryu, 
Jones, and Frank, 2000) is involved. See for example, Diver, Brown, and 
Rust (1996) who proposed that the KHI evolves into centers for the tearing 
mode instability, which would account for the Ellerman bombs that 
frequently appeared above the MDFs. 

It is interesting to contrast the properties of the inflowing MDFs with 
the dipolar outflowing MMFs previously reported around decaying 
sunspots. Outflowing MMFs move faster than the surrounding granules; 
they may appear at first to be a continuation of penumbral fibrils or they 
may develop anywhere in the moat; the two poles separate gradually; and 
they disappear when they reach the sunspot moat, or they just fade away 
(Shine and Title, 2001). The pole nearest the spot usually has the same 
polarity as the spot (Harvey and Harvey, 1973). However, Yurchyshyn, 
Wang, and Goode (2001) observed some outflowing dipolar MMFs with 
the nearest pole polarity opposite to that of the parent spot. This is the 
same orientation that we found for die inflowing MDFs. The inflowing 
MDFs that we observed moved with the surrounding material; they were 
associated with the dark fibrils of the emerging flux; they could develop 
anywhere in the emerging flux region. When an MDF reached the 
negative-polarity sunspot, the positive part dissolved and the spot 
developed a peninsula of negative field that stretched out to incorporate 
the negative part, which trailed the positive one. 

Ryutova et al. (1998) and Harvey and Harvey (1973) suggested that 
outflowing dipolar MMFs are convex kinks in the horizontal fields 
extending outside penumbrae. From their shape, they are called Q-loops. 
The sunspot-decay MMFs in our own observations would appear to fit this 
description. Near the northern edge of the trailing spot we saw one MMF 
flowing outwards, and indeed, its vector field configuration revealed an Q- 
loop. However, there seems no escaping the conclusion that the sunspot- 
building MDFs are U-loops formed between segments of emerged flux. 
The MDFs cannot be kinks in the rising flux ropes because then we would 
see the parts with opposite polarity separating. Instead, the two MDF 
parts formed by pulling apart, but then they remained connected for the 
entire life span of the MDF. We conclude that they formed when plasma 
on just emerged field lines cooled, then pooled, then distorted the field 
lines, making a concave depression in them. The distorted segment, a U- 
loop, is swept by the convective flows to the edge of a sunspot or 
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supergranule. There the entrained material can slide down one of the main 
legs of the flux rope bundle in the active region, which presumably 
extends deeply into the convection zone, possibly all the way to the base. 

These observations illustrate how knots of plasma entrained in 
segments of emerging flux ropes are propelled along the solar surface, like 
pearls on a string, towards sunspots and the chutes at supergranules 
boundaries. There, they can slide back down into the sea of plasma from 
which they came. The observations also illustrate dramatically how vector 
magnetograms facilitate the interpretation of active region phenomena. 

Our results may have an application that is wider than active regions - 
there is a background pattern of dipolar fields everywhere on the sun. 
Dipolar features are often bom with the same pattern as those seen here, 
namely, the positive and negative poles are bom together, then separate to 
assume a characteristic distance from each other. Since our vector 
magnetic field observations showed that the horizontal field in the MDFs 
is far more important than the vertical component, we suggest that there 
could be a background pattern of horizontal fields everywhere. They 
could reside in the convectively stable layer just above the granulation. 
The plasma sloshing around in those fields should often form pools and 
deform those fields, giving the appearance in longitudinal magnetograms 
of magnetic dipoles coming and going. The difficult-to-understand 
background fields (Zhang et al., 1998) may, in part, be due to dips in 
horizontal fields not observed, so far. Not until there are very sensitive 
vector field measurements of the weak fields outside active regions will 
we know the true structure of the ephemeral dipolar features. 
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ABSTRACT 

We investigate the statistical properties of Ellerman bombs in the dynamic emerging flux 
region NOAA 8844, underneath an expanding arch filament system. High-resolution chromo- 
spheric Ha filtergrams (spatial resolution 0.8 ,/ ), as well as photospheric vector magnetograms 
(spatial resolution 0.5"), and Dopplergrams, have been acquired by the balloon-borne Flare Gen- 
esis Experiment. Ha observations reveal the first “seeing-free” data set on Ellerman bombs and 
one of the largest samples of these events. We find that Ellerman bombs occur and recur in 
preferential locations in the low chromosphere, either above, or in the absence of, photospheric 
neutral magnetic lines. Ellerman bombs are associated with photospheric downflows and their 
loci follow the transverse mass flows on the photosphere. They are small-scale events, with typi- 
cal size 1.8 x 1.1 arcsec ? , but this size depends on the instrumental resolution. A large number 
of Ellerman bombs are probably undetected, due to limited spatial resolution. Ellerman bombs 
occur in clusters which exhibit fractal properties. The fractal dimension, with an average value 
- 1.4, does not change significantly in the course of time. Typical parameters of Ellerman bombs 
are interrelated and obey power-law distribution functions, as in case of flaring and sub-flaring 
activity. We find that Ellerman bombs may occur on separatrix, or quasi-separatrix, layers, in the 
low chromosphere. A plausible triggering mechanism of Ellerman bombs is stochastic magnetic 
reconnection caused by the turbulent evolution of the low-lying magnetic fields and the contin- 
uous re-shaping of separatrix layers. The total energies of Ellerman bombs are estimated in the 
range (10 ,10 erg , the temperature enhancement in the radiating volume is 2 x 10 1 * 3 K and 
the timescale of radiative cooling is short, of the order of a few seconds. The distribution function 
of the energies of Ellerman bombs exhibits a power-law shape with an index ~ -2.1. This sug- 
gests that Ellerman bombs may contribute significantly to the heating of the low chromosphere 
in emerging flux regions. 

Subject headings: Sun: chromosphere — Sun: flares — Sun: magnetic fields — Sim: photosphere — Sun: UV 
radiation 


1. Introduction 

Ellerman (1917), in a pioneering work, reported 

a solar phenomenon “so extraordinary that it 
seemed hardly real”. Spectroscopic studies at the 
Mount Wilson Solar Observatory had discovered a 
brilliant emission kernel at both wings of the Ha 


line in several locations of the low chromosphere. 
These short-lived, small-scale events are now usu- 
ally called Ellerman bombs (EBs). 

The spectral profile of EBs consists of a deep 
absorption core at the Ha line center with ex- 
cess emission at the red and the blue wings of the 
line. The emission reaches a maximum at wave- 
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lengths ~ Ha =t 1 A and fades at typically (N/ ±5 
A from the line center (Engvold & Maltby 1968; 
Severny 1968). An asymmetry with excess emis- 
sion at the blue wing has been first reported by 
Koval & Severny (1970) and Bruzek (1972). Be- 
cause of the peculiar appearance of bomb spectra 
on photographic negatives, EBs are also known 
as “moustaches”. The first to confirm the iden- 
tity between EBs and moustaches were McMath, 
Mohler, & Dodson (1960). 

For decades, EBs have been known as frequent, 
short-lived brightenings, seen mostly at the wings 
of the Ha line. They generally occur in areas of 
flux emergence and strong magnetic fields, such 
as the vicinity of the penumbra around isolated 
sunspots, and the interspot areas of new and dy- 
namic emerging flux regions, beneath the expand- 
ing arch filament system (AFS). Zachariadis, Alis- 
sandrakis, & Banos (1987) estimated a mean oc- 
currence rate of about 1.5 EBs min -1 in an area of 
size 18 x 24 arcsec 2 . The typical sizes of EBs are 
of the order one arcsec. Mean lengths of 1 - 1.4 
arcsec have been reported. In most cases, EBs 
appear as elongated structures with length typi- 
cally 37 - 45% larger than their width (Kurokawa 
et al. 1982; Nindos & Zirin 1998). The typical 
EB lifetimes were estimated from early works to 
be ~ 20 min (Severny 1965). Bruzek (1972) and 
Vorpahl and Pope (1972), however, reported that 
certain EBs could be visible for times exceeding 
1 hr. The estimation of EB durations depends 
critically on the temporal resolution (cadence) of 
the observing instrument, since EBs often tend to 
recur at approximately the same locations (Vor- 
pahl & Pope 1972; Bruzek 1972; Zachariadis, Alis- 
sandrakis, & Banos 1987). The generally accepted 
mean lifetime of EBs is 10 - 14 min (Roy & Lep- 
arskas 1973; Kurokawa et al. 1982; Payne 1993; 
Nindos & Zirin 1998). In a recent, high-cadence 
(~ 1 s) series of observations, Qiu et al. (2000) 
reported mean EB lifetimes of 10 - 20 min. 

In studies of the surroundings of EBs, Severny 
(1968) detected intense ascending chromospheric 
motions and Rust (1968) found that EBs occur 
above “satellite” (opposite polarity) sunspots seen 
at the vicinity of well-defined spots. EBs are 
also associated with Ha surges, whose threads 
(~ 5 x 10 3 km) follow the magnetic field lines 
calculated from current-free extrapolations (Rust 
1968). Several EBs were found at the boundaries 


of well-defined magnetic structures (Roy and Xep- 
arskas 1973; Dara et al. 1997), related to dipo- 
lar emerging flux (Zachariadis, Alissandrakis, & 
Banos 1987), and moving magnetic features (INin- 
dos and Zirin 1998). 

Despite the fact that EBs have long been 
known, their triggering mechanism remains puz- 
zling. In early studies by Severny and co- workers 
(for a review see Severny 1964), EBs, or mous- 
taches, were viewed as low-altitude sub-flares. 
This is still a highly debated viewpoint, with sev- 
eral issues yet to be resolved. One of these issues is 
the temporal profiles of EBs. In contrast to dares 
and sub-flares, the temporal profiles of EBs are 
rather symmetric (Payne 1993; Qiu et al. 2000). 
Flares and sub-flares, on the other hand, mostly 
show an impulsive, short rise phase, followed by 
a much longer, gradual decay phase. Other EB 
properties include association with facular gran- 
ules (Bruzek 1972; Dara et al. 1987), Ha filigree, 
calcium bright points and white-light faculae ("Wil- 
son 1981), facular points (Kitai & Muller 1984), 
photospheric filigree (Denker 1997), and moving 
magnetic features (Nindos and Zirin 1998). Vor- 
pahl & Pope (1972), found that EBs correlate 
with photospheric bright points (BPs) at 3840 A , 
which reach maxima at ~ 3.5 ±1.0 min prior to 
EB occurrence. Qiu et al. (2000) showed that EBs 
can also be seen as BPs in the Ultraviolet CUV) 
continuum emission at 1600 A. These events occur 
in the temperature-minimum (T min ) region, at a 
height difference ~ 100 km from Ha EBs. 

Several models have been proposed to accom- 
modate both the morphology of EBs and “their 
spectral features. None of these models is gen- 
erally accepted. Severny (1968) proposed a model 
with an expanding opaque layer that exhibits a ve- 
locity gradient with height, in a turbulent medium 
characterized by an ambient macro-turbulenrt ve- 
locity. A suitable choice of parameters could ac- 
count for the distribution of intensity wingjs in 
moustaches, the blue asymmetry and the absorp- 
tion gap at the Ha line center. Pikel’ner (L.974) 
argued that if EBs are triggered by magnetic re- 
connection, then the temperature should lo-cally 
reach coronal values after field-line annihilation. 
Kitai (1983) found that a hot (AT ~ 1.5 x lO 3 K) 
and dense (^ ~ 5) gradient in the atmosphere 
may account for the line broadening in EBs. IKitai 
& Muller (1984) viewed EBs as events occu iring 
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in funnel-like, individual flux tubes, rather than as 
results of magnetic reconnection. Stellmacher & 
Wiehr (1991) attributed the moustache effect to 
the superposition of a hot component represent- 
ing the magnetic structure, surrounded by a cool 
intergranular region. Hu, Song, & Li (1995) con- 
structed the first magnetohydrodynamic (MHD) 
model of EBs. The model was based on mag- 
netic reconnection on a current sheet, formed by 
the motion of an ascending plasmoid triggered by 
ponderomotive force. The local temperature in 
EBs reached 10 4 K. Higher temperatures were 
avoided, since the field-line annihilation products 
were carried away instantly by suprathermal par- 
ticles. A similar MHD model, modified for a 
plasma with low ionization degree, was proposed 
by Ji, Song, & Li (2001). Ding, Henoux, & Fang 
(1998) and Henoux, Fang, & Ding (1998) pro- 
posed a non-thermal EB model that attributed 
the moustache line profiles to the bombardment 
of the low chromosphere by energetic particles. 
These particles were expected to originate from re- 
connection sites higher in the atmosphere. Diver, 
Brown, & Rust (1996) attributed EBs to a tearing- 
mode instability which was initially triggered by 
a purely hydrodynamical Kelvin-Helmholtz insta- 
bility. The dissipated Ohmic power was estimated 
at ~ 10 21 er^ s"" 1 . Other energy estimates of EBs 
are in the range (10 25 , 10 27 ) erg (Bruzek 1972; Hu, 
Song, & Li 1995; Henoux, Fang, & Ding 1998). 

To address the outstanding problem of the trig- 
gering cause of EBs, we performed a detailed 
analysis of EB emergence, morphology, statis- 
tical properties, and energetics. We used the 
first “seeing-free” sample of EBs, obtained by the 
balloon-borne Flare Genesis Experiment (FGE - 
Rust 1994a; Bernasconi et al. 2000; Bernasconi, 
Rust, & Eaton 2001). Several hundreds of EBs 
appear in our data, therefore making our sample 
adequate for statistical treatment. 

The paper proceeds as follows: In § 2 we dis- 
cuss the acquisition of the FGE data and the data 
reduction process. In § 3 we estimate the mean 
occurrence rate of EBs and describe their spatial 
distribution. In § 4 we discuss the morphology of 
EBs and their relation to the underlying photo- 
spheric conditions and to activity in the T min re- 
gion. In § 5 we perform the first statistical study 
of EBs, based on the fractal properties of the EB 
clustering process and on typical EB parameters, 


such as the spatial extent, peak and total activity. 
In § 6 we estimate EB energies relying on estima- 
tions of radiative losses and mass motions in lEBs. 
In § 7 we propose an interpretation of EBs, sup- 
ported by our observations. In § 8 we summarize 
our findings and discuss the implications of our 
results. 

2. Flare Genesis Experiment and Active 
Region NOAA 8844 

The data used in this work were collected dur- 
ing the second FGE flight, in January 2000. The 
flight duration was 17 days (January 10 to Jaituary 
27). The apparatus was installed on a balloon- 
borne gondola and it circumnavigated the Antarc- 
tic continent at an altitude of 37 km , obtaining 
an un-obstructed, 24 - h view of the solar disk 1 . 
The principal optic of the FGE was an 80 an 
Cassegrain telescope (the largest solar telescope 
ever flown) with a F/ 1.5 light-weight, ultrat-low 
expansion primary mirror, a single-crystal silicon 
secondary and a graphite-epoxy telescope tube. 
The tube structure provided high thermal stabil- 
ity and the silicon secondary mirror gave optimum 
heat transfer with minimum thermal distortion. 
For this reason, we believe that the acquired im- 
ages are free from distortion due to turbulence in 
the telescope, or in the atmosphere (Bernasconi et 
al. 2000; Bernasconi, Rust, & Eaton 2001). The 
focal plane instrumentation consisted of a polar- 
ization analyzer unit (two liquid-crystal polariza- 
tion modulators and a linear polarizing filter), a 
selection of 1.25 A order-isolation filters, a 0.16 
A tunable Fabry-Perot etalon filter (Rust 1994b) 
and refractive optics to record a 160 — arcsGC Q) 
section of the Sun on a 1024 x 1024-pixel OCD 
camera. The field of view of the CCD was 92 x 
92 arcsec 2 , the plate scale was 0.091 arcsec. per 
pixel, and the achieved image quality (spatial res- 
olution) varied between 0.5 arcsec and 0.8 ar csec , 
depending on the observed wavelength. To achieve 
a better signal-to-noise ratio, the images werre in- 
terpolated into a uniform grid of 512 x 512 pixels, 
with pixel size equal to 0.18 arcsec . 

FGE collected nearly simultaneous observations 
of the photospheric magnetic field vector, the line- 

*For a complete description of the FGE missiozn see 
Bernasconi et al. (2000), as well as the FGE we"b site 
http://sd-www.jhuapl.edu/FlareGenesis / 


3 



of-sight (LOS) longitudinal photospheric velocity 
(from Dopplergrams) and filtergrams at the blue 
wing of the Ha line. Let us discuss briefly each of 
these data sub-sets: 

(1) The Stokes parameters required for vector 
magnetograms were obtained from observa- 
tions at the red and the blue wings of the 
photospheric Ca I line (6122.2 A). They 
were obtained from sequentially recorded 
sets of 12 filtergrams at the polarization 
states I + Q, I -Q, I + U, I -U, I + V 
and I — V at ±0.08 A at the wings of the 
Ca I line. The Stokes images (Q/I), ( U/I ) 
and (V/I) were calculated from the normal- 
ized difference between measurements of op- 
posite polarization states. The magnitudes 
of the LOS longitudinal and transverse mag- 
netic field components were calculated fol- 
lowing the method of Ronan, Mickey, & Or- 
rall (1987). The precise location of the AR 
on the solar disk was found by comparison 
to simultaneous, full-disk Kitt Peak magne- 
tograms. This enabled the calculation of the 
local reference system (LRS) photospheric 
magnetic field components. The resolution 
of the 180°-ambiguity was qualitative, rely- 
ing on a combination of a potential field solu- 
tion and an empirical, expected orientation 
of the transverse component (Bemasconi et 
al. 2002). 

(2) Dopplergrams were obtained by calculat- 
ing the normalized difference between filter- 
grams obtained at the wings of the Ca I line. 
During each 12-image set, an averaging pro- 
cess was performed. Standard filters were 
applied to remove some effects of the 5 min r 
oscillation and to cancel out image-to-image 
intensity variations introduced by vector po- 
larimetric measurements. 

(3) Filtergrams at the blue wing of the Ha line 
(6562.8 A) were obtained between each set 
of Stokes images. Corrections for the dark 
current and the gain function were applied 
and the filtergrams were carefully co-aligned 
by means of standard cross-correlation tech- 
niques. During the co-alignment process, 
the center of the AR was kept fixed at the 
center of the field of view. 


Restrictions imposed by the observing plan, the 
instrument, the CCD detector, and the data stor- 
age medium resulted in a cadence varying in the 
range (3.5, 7.5) min. The Ha data refer solely to 
the blue wing of the line, at about 0.8 A from the 
line center. For these reasons, we will not attempt 
any conclusions regarding either the temporal or 
the spectral profiles of EBs. 

On January 25, 2000, FGE was observing the 
young, dynamic, emerging flux region NOAA 
8844, just after the active region (AR) shewed 
signs of rapid development. An expanding AFS 
was observed in Ha and a multitude of EBs was 
detected beneath the large-scale structure. Axi un- 
interrupted data set was obtained for 3.5 h, (UT 
15:50 - 19:16). The AR entered its declining phase 
24 h after the observations. The AR’s mean heli- 
ographic coordinates were 5° N 30° W. The axis 
of the AR formed a tilt single of /v 28° with the 
equator. The AR consisted of a compact pstir of 
small sunspots which were moving apart rapidly. 
Although the AR produced no flares, some sub- 
flaring activity was reported by the NOAA Space 
Environment Center, with the most prominent 
sub-flares observed at 16:02, 16:27, 17:04 and 18:54 
UT. 

In this work, we will use the following FGE 
data: 

(i) 55 vector magnetograms with a spatial 
resolution ~ 0.5 arcsec and a cadence 
~ 3.5 min. 

(ii) 28 Dopplergrams with a spatial resolution 
~ 0.8 arcsec and a cadence ~ 7.5 min , 

(iii) 28 Ha — 0.8 A filtergrams with a spa- 
tial resolution ~ 0.8 arcsec and a cadence 
~ 7.5 min. 

The FGE data on AR NOAA 8844 were 
complemented by simultaneous observations by 
Yohkoh/ SXT, SOHO/MDI, TRACE and the 
Imaging Vector Magnetograph (IVM) of the Uni- 
versity of Hawaii. We have not yet fully retrieved 
this multi-instrument data collection. In. this 
study, however, we make a first step in coupling 
the FGE data with a data set from TRACE which 
refers to UV 1600 A continuum emission im the 
'min region. We have selected 13 TRACE frames 
that are simultaneous with respective FGE frames. 
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Fig. 1.— Four nearly simultaneous views of the AR NOAA 8844 on January 25, 2000, obtained by the FGE: 
(a) LOS longitudinal photospheric magnetogram. Bright /dark patches indicate positive/negative polarity. 
The white ellipse encloses a supergranular cell, at the vicinity of the following sunspot, (b) Photospheric 
Doppler velocities. Bright/dark patches indicate upflows/downflows. c) tfa-0.8A filtergram. A multitude 
of EBs may be seen as localized brightenings underneath the superposed AFS (elongated, dark fibrils) (d) 
UV 1600 A continuum emission, obtained by TRACE. The tick mark separation in all images is 1 arcsec 
North is up; west is on the right. 



The co-alignment of the TRACE 1600 A images 
with the FGE data was accomplished in two 
stages: First, by aligning the position of the 
sunspots in both wavelengths and, second, by 
aligning small-scale features, visible in both wave- 
lengths. To map the transverse velocity field on 
the photospheric level (proper motions), we have 
applied a cross-correlation tracking technique to 
FGE white light images (R. A. Shine, 2001, pri- 
vate communication). 

A set of nearly simultaneous views of AR 
NO A A 8844 is shown in Fig. 1. A typical LOS 
photospheric magnetogram is shown in Fig. la. 
The photospheric Doppler velocity map in shown 
in Fig. lb. The simultaneous Ha - 0.8 A fil- 
tergram is given in Fig. lc. Notice the numerous 
EBs beneath the AFS (dark elongated structures 
extending between the two sunspots). EBs are 
visible as localized enhancements of Ha — 0.8 A 
emission. In Fig. Id, the simultaneous 1600 A 
continuum image of the T min region is given. Lo- 
calized enhancements of emission are also abun- 
dant here. We will, hereafter, refer to these tran- 
sients as 1600 A BPs. Besides BPs in Fig. Id, 
we also obtain bright, elongated structures in the 
interspot area. We believe that these fibril struc- 
tures originate from emission in the C IV line 
(1550 A). This emission is beyond the scope of 
the present study. 

Figure 1 shows that the emergence of mag- 
netic flux in the AR is an asymmetric process. 
The preceding sunspot (positive polarity) is more 
compact than the trailing sunspot (negative po- 
larity), which appears less organized and more 
fragmented. Westward motions at the vicinity of 
the leading spot are more systematic compared to 
those at the vicinity of the following spot, where 
we notice a small (~ 15 Mm) supergranular cell 
(enclosed by the ellipse in Fig. la). We will not 
deal with these asymmetries in this work, but one 
may refer to the works by van Driel-Gesztelyi and 
Petrovay (1990), Petrovay et al. (1990), and Fan, 
Fisher, & DeLuca (1993). 

3. Occurrence Rates and Spatial Distribu- 
tion of EBs 

To identify individual EBs, we rely on the con- 
trast J c (x,£) = , calculated for each loca- 

tion in the Ha — 0.8 A filtergrams. Here /(x,t) is 


the intensity of a pixel with vector position x at 
time t and Iq is the background (mean) intensity 
for all frames and for the entire observing time. 
Setting all pixels with contrast fainter than some 
predefined value to zero, we obtain “islands” of 
enhanced Ha — 0.8 A emission, which correspond 
to EBs. 

To track the temporal evolution of EBs we as- 
sume that an EB, identified in one frame, survives 
to the next frame if there is a spatial overlapping 
of at least one pixel between islands of enhanced 
Ha — 0.8 A emission in two consecutive frames. 
This somewhat crude rule may underestimate the 
total number of EBs and overestimate EB lifetimes 
because: (i) Experiments with better cadence re- 
port EB lifetimes of the order of a few min, which 
is less than our cadence. The tendency of EBs to 
recur at the same locations may, therefore, lead 
to undetected EB recurrence. ( ii) Depending on 
the local proper motions, an EB may migrate to 
the location of a pre-existing EB which may have 
faded between frames. We have created an algo- 
rithm that automatically identifies EBs and tracks 
their evolution in time. According to the algorith- 
mic rules, (i) two EBs that merge are, thereaifter, 
treated as one event, and (ii) an EB that splits 
in two or more separate emission kernels is still 
treated as one EB and is considered terminated 
only when the last remaining kernel fades t>elow 
background. 

The results of the above selection process are 
the following: 

The observed contrasts of EBs vary in the range 
[0.05,0.30], that is, 5% to 30% above the Back- 
ground. The average contrast, calculated for all 
pixels over the entire observing time, is ~ 0.08. 
This average is significantly lower compared to 
previous reports. Zachariadis, Alissandrakis, & 
Banos (1987) reported an average contrast of ~ 
27% , while Nindos & Zirin (1998) reported a mean 
value of ~ 29 — 31%. We have defined several 
contrast thresholds, or background contrasts, in 
the range [0.05,0.20]. The total number of EBs 
varies between 593 events (for a background con- 
trast I Co = 0.05) and 29 events (for I Co — 0.20). 
For I Co = 0.08, close to the average contrast, 
we identified 295 EBs. The total observing time 
being ~ 206 min and assuming a smooths un- 
interrupted EB production, we find that the mean 
occurrence rate of EBs varies in the range 
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Fig. 2. The mean origin locations of EBs (crosses) obtained using a background contrast I Co = 0.08 
(295 events). A typical LOS longitudinal magnetogram is shown in grayscale. Bright/dark patches indicate 
positive/negative polarity. The EB-producing area has been divided in three parts (for details, see text). 
The three arrows indicate the locations of the EB examples shown in Fig. 5, one for each of the designated 
areas 1, 2, and 3. The tick mark separation is 1 arcsec. North is up; west is on the right. 


[0.14,2.88] EBs min 1 , for an EB producing 
area ~ 1800 arcsec?. For I Co = 0.08, the mean 
occurrence rate is ~ 1.43 EBs min -1 . The above 
occurrence rates for a given background contrast 
should be considered as lower limits of the actual 
occurrence rates, since (i) our results may be sys- 
tematically biased due to the moderate cadence, 
compared to typical EB lifetimes, and (ii) EBs oc- 
curring directly underneath the AFS’s dark fibrils 
are obscured by the large-scale structure. 

In Fig. 2 we show the mean locations of origin 
for the 295 EBs detected for I c o ~ 0.08. A typical 
LOS longitudinal magnetogram is the grayscale 
background. It is clear from Fig. 2 that EB pro- 
duction is inhomogeneous and enhanced at: 


(1) The boundaries of the supergranular cell, at 
the vicinity of the following sunspot (area 
1). The supergranular boundaries are cov- 
ered by EBs, while EB density is signifi- 
cantly lower in the interior of the super- 
granule. The boundaries of the cell are 
sustained by converging proper motions, 
of the order (0.2, 0.5) km s -1 , and they 
are dominated by downflows in the range 
( — 0-4, —0.1) km a -1 . In the supergxanu- 
lar interior, EBs are mostly associated with 
small-scale moving dipolar (magnetic) fea- 
tures (MDFs) which drift toward the follow- 
ing sunspot, or toward the cell boundaries, 
with velocities in the range (0.7,1) km s -1 . 
MDFs found in the FGE data are peculiar 
moving magnetic features and are studied 
separately by Bemasconi et al. (2002) . 
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(2) The polarity reversal area (area 2). EBs ob- 
served in this area of intense flux emergence 
are often associated with small-scale MDFs 
(Bernasconi et al. 2002) moving with veloc- 
ities of the order (0.3, 0.8) km s -1 . 

(3) The vicinity of the leading sunspot (area 
3). This area is mostly unipolar and char- 
acterized by systematic proper motions of 
newly emerged magnetic fields of positive 
polarity toward the sunspot. Proper mo- 
tions are of the order (0.6, 0.8) km s~ x and 
the moving patches correlate with downflows 
~ —0.5 km s" 1 . Notice that the origin lo- 
cations of EBs form chains that lead to the 
preceding sunspot. These chains correlate 
well with the paths of the above-mentioned 
moving bundles of magnetic flux. 

Our results are in agreement with previous studies 
reporting that the spatial distribution of EBs is 
not random or homogeneous (Vorpahl and Pope 
1972; Bruzek 1972; Zachariadis, Alissandrakis, & 
Banos 1987). A detailed study of the photospheric 
patterns for each of the three areas shown in Fig. 
2, follows in § 4.2. 

4. Morphology and Characteristics of EBs 

4.1. Typical size and the degree of elonga- 
tion 

Previous works (Kurokawa et al. 1982; Nindos 
and Zirin 1998) reported that EBs are elongated 
structures with their length typically 37% - 45% 
larger than their width. To investigate this, we 
selected five representative background contrasts 
J Co , namely, 0.05, 0.08, 0.12, 0.16, and 0.20 and 
we inspected visually the shapes of the most con- 
spicuous EBs for each case. We measured the sizes 
of 328 EBs (165 for I Co = 0.05, 103 for I Co = 0.08, 
40 for I C0 = 0.12, 15 for I CQ = 0.16, and 5 for 
I co = 0.20). The main/secondary axis of the ex- 
amined EBs was called the length/width of the 
EBs. Knowledge of the length l and the width w 
of an EB allows the calculation of the degree of 
elongation e, defined as e = x 100. From this 
definition, e = 0% for a circle and e ~ 100% for 
l > w. 

In Fig. 3 we show the average degree of elon- 
gation for all five definitions of the background. 



0 10 20 30 40 50 60 70 SO 
Elongation (%) 


Fig. 3. — Histogram of the measured degree 

of elongation for well-defined EBs obtained with 
background contrasts Ic 0 = 0.05, 0.08, 0.12, 0.16 
and 0.20. Also shown, with error bars, is the stan- 
dard deviation for each EB percentage associated 
with a certain degree of elongation. 

The histogram was created using 8 equal elon- 
gation bins of width 10%. The measured degree 
of elongation is located at the center of eacki bin 
while the error bars correspond to the stam dard 
deviation for each EB percentage associated with 
a degree of elongation. From Fig. 3 we con- 
clude that (i) EBs are elongated structures, but 
there is an upper limit for the degree of elongation 
and (it) there exists a preferential degree of elon- 
gation. Roughly 10% of EBs are circular struc- 
tures (e ~ 0), while for 75% of EBs we obtain 
10% < e < 40%. For nearly half of the elongated 
EBs we obtain e ~ 40%, that is, their length is 
40% larger than their width. This is the prefer- 
ential degree of elongation, even given the error 
bars. This degree of elongation is in close agree- 
ment with previous results. For about 15% of EBs 
we obtain e > 40%, while no EB exhibits a degree 
of elongation larger than ~ 70%. We conclude 
that EBs are mostly elongated structures, but only 
rarely (if at all) does the degree of elongation ex- 
ceed 60%. 

The mean EB sizes, on the other hand, can be 
estimated by more than one way: 

(a) By directly measuring EB sizes. This 
method is biased in favor of the most con- 
spicuous EBs, since the measurements can- 
not be reliable for EBs with sizes close t o the 
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spatial resolution. Small EBs, consequently, 
have to be excluded and the mean EB sizes 
are calculated from only a fraction (less than 
50%) of the total event sample. 

(b) By calculating EB sizes. Using the prefer- 
ential degree of elongation, we assume that 
l ~ 1.67m and that EBs are elliptic al struc- 
tures. We then obtain w ~ where S 

is the measured EB area. This set (/, w) of 
EB dimensions can be viewed as the antic- 
ipated dimensions of an EB which extends 
over an area S. 

We noted in § 3 that the average contrast of 
our EBs is ~ 0.08. Let us use this background 
(J Co = 0.08) in order to determine the mean EB 
sizes. By directly measuring the EB sizes, that 
is, by using 103 events out of a total of 295 
EBs, we obtain average EB sizes of the order 
2.5 x 1.7 arcsec 2 . Assuming l ~ 1.67 w, we use 
the entire EB sample to obtain an average size of 
the order 1.8 x 1.1 arcsec 2 , using the measured 
EB areas. The mean EB dimensions are not much 
larger than our spatial resolution (~ 0.8 arcsec). 
Improved spatial resolution will, probably, result 
in even smaller “typical” EB sizes. This is an in- 
dication that finite spatial resolution precludes a 
detailed analysis on EB sizes and that probably 
many EBs remain undetected. This assertion is 
further supported by the significantly larger av- 
erage contrast in Ha - 0.8 A within the EB- 
producing area (Fig. lc), excluding EBs, com- 
pared to the average contrast of the quiescent ar- 
eas on both sides of the interspot region. More- 
over, we will show in §5.2 that there is a strong 
positive correlation between the size (area) and 
the peak intensity of EBs, indicating clearly that 
fainter EBs occupy smaller areas. 

4.2. The photospheric signature of EBs 

The photospheric conditions that give rise to 
EB occurrence are less than clear. Topics to be 
addressed include the association of EBs with neu- 
tral lines of the longitudinal photospheric mag- 
netic field, photospheric Doppler motions, and the 
proper motions of EBs. In the following analysis 
we investigate the correlations of EBs with the un- 
derlying photospheric pattern. We will, again, use 
a background contrast I Co = 0.08, close to our av- 


erage EB contrast. 

We first study the association of EBs with pho- 
tospheric neutral lines of the LOS longitu<iinal 
magnetic field. In this analysis we do not use 
the tracking algorithm of EBs in time, since the 
process would then become very complicated.. In- 
stead, we examine each Ha filtergram, counting 
each EB as a different case. Consequently, we 
examined 173 EBs in area 1, 172 EBs in area 2, 
and 340 EBs in area 3. Our results suggest that 
the presence of photospheric neutral lines is not 
a prerequisite for EB triggering. Area 1 is domi- 
nated by persistent neutral lines that outline the 
supergranular boundary. Since these boundaries 
are areas of enhanced EB production, a decisive 
majority of EBs (77%) is found to occur above 
neutral lines. This percentage includes EBs that 
occur above the neutral lines of MDFs in th.e su- 
pergranular interior. Area 2 is the polarity rever- 
sal area of the AR. The majority of EBs found 
in this area (65%) also occur above neutral lines 
of the newly emerged magnetic flux. These lines 
are not as well-defined as the neutral lines in area 
1, because magnetic fields are weaker in area 2. 
In some cases, we see EBs occurring above small- 
scale MDFs. Area 3, on the other hand, is mostly 
unipolar. Nonetheless, this area also produces a 
large number of EBs. Only 36% of these EBs cor- 
relate with photospheric neutral lines, while the 
majority of EBs (64%) occur in unipolar areas. 

Next, we investigate the proper motions of EB 
loci. To measure these horizontal velocities, we 
calculate the mean locations (centroids) of EBs for 
the first and the last frame in which they appear. 
The mean velocity of an EB loci, then, is the ratio 
of the distance covered by the EB’s centroid, over 
the elapsed time between the initial and the final 
frame in which the EB was detected. 

Our results suggest that the supergrauxiular 
boundaries (area 1) trigger numerous EBs that 
are mostly short-lived and their loci do not travel 
significant distances. The calculated proper mo- 
tions do not exceed 0.4 km s -1 , while, in most of 
cases, they are less than 0.2 km s -1 . In the snper- 
granular interior the situation is different, with a 
few EBs (associated with MDFs) showing proper 
motions in the range (0.4, 0.8) km s -1 , or even 
exceeding 0.8 km s~ 1 . The polarity reversal area 
(area 2) produces EBs with typical loci proper 
motions in the range (0.2, 0.4) km s -1 . 
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Fig. 4. Two examples of EB proper motions in area 1 (white rectangle) and in area 3 (black rectangle). 
The evolution of EBs in three consecutive frames, labeled (a), (b), and (c), is shown. The un-scaled image 
is a typical LOS longitudinal magnetogram. The enlarged areas include the LOS longitudinal magnetic 
field in grayscale, the EB contrast (contours), and the average proper motions (vector field), calculated by 
correlation-tracking techniques. Tick mark separation in the magnification is 1 arcsec. The upper area 
(white rectangle) includes a strong neutral line, while the lower area (black rectangle) is mostly unipolar. 
The EB locus in area 1 moves with a velocity ~ 1.15 km a -1 , whereas the EB locus in area 3 moves with 
velocity ~ 1.35 km s 1 . The UT times are shown at the left comer of the enlarged images. 
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Slow-moving EBs (with velocities < 0.2 km s -1 ) 
were also found. In only one case do EB proper 
motions exceed 0.4 km s -1 in area 2. The vicinity 
of the leading polarity (area 3), on the other hand, 
produces a large number of EBs with loci showing 
significant proper motions. The typical velocity 
range is (0.4, 0.6) km s -1 . Fast motions, of the 
order (0.6, 0.8) km a -1 , were found in some cases. 
Several EBs move with velocities < 0.2 km s -1 , 
while very few EBs show velocities that exceed 
0.8 km a -1 . 

EB loci appear to follow the transverse photo- 
spheric mass flows. This is more definitive in the 
areas of systematic horizontal flows, i.e., the su- 
pergranular area and the unipolar vicinity of the 
leading sunspot. We depict two examples in Fig. 
4. Two parts of the field of view have been en- 
larged, in area 1 and in area 3, respectively. In the 
expanded images, we depict the LOS longitudinal 
magnetic field (grayscale), the EBs (contours) and 
the average proper motions (vector field). Three 
frames, labeled (a), (b), and (c), show the tempo- 
ral evolution of EBs, in both examples. Clearly, 
the EBs, or at least the EB lod, appear to follow 
the photospheric flow pattern. 

Let us now study the association of EBs with 
the LOS longitudinal velocity of the photospheric 
magnetic fields (Fig. lb). Before describing 
the results, it is useful to give the overall pic- 
ture of the Doppler motions in the AR. Flux 
emergence is asymmetric, as noted in § 2. The 
sunspots are mostly associated with upflows. In 
the leading spot, ascending motions span in the 
range (0.1, 0.4) km s -1 , whereas in the follow- 
ing spot we obtain intense upflows of the order 
(0.2, 1.0) km a -1 . No preferential longitudinal 
motions of magnetized plasma are detected in the 
polarity reversal area. Doppler velocities vary in 
the range (—0.3, 0.3) km s -1 . Persisting down- 
flows, on the other hand, can be found (i) at the 
supergranular boundaries, with motions of the or- 
der (-0.4, —0.1) km a -1 , and (ii) at the loca- 
tions of the moving unipolar magnetic patches, 
in area 3. Doppler velocities in this area are 
~ —0.5 km s~ l . 

To associate EBs with the photospheric Doppler 
motions, we performed a similar analysis as in the 
case of EB association with neutral magnetic lines. 
We find that almost all EBs are associated with 
photospheric downflows. In area 1, we find that 


photospheric downflows accompany 86% of EBs. 
The respective percentages are 86% for area 2 
and 81% for area 3. Notice the similar percent- 
ages in all areas of the AR. This behavior is not 
seen either in the association of EBs with neutral 
lines or in the proper motions of EB loci. Dif- 
ferent areas of the AR give different results, in 
both cases. We need to clarify that, to find the 
percentage of EBs associated with downflows in 
area 1, some further considerations were made. 
For several EBs in area 1 the Doppler signal was 
blue-shifted. However, in most of these cases, the 
Doppler velocity exhibited a local minimum in ab- 
solute magnitude (examples are given in § 4.4). 
These relative downflows in the locations of EBs, 
with respect to the EBs’ surroundings, may be 
either due to LOS viewing angle problems, or due 
to small-scale, unresolved downflows superposed 
to macroscopic upflows. In area 1, 46% of EBs 
are associated with downflows and 40% of EBs 
correlate with upflows showing a clear velocity 
minimum in the EB location. Adopting the above 
interpretation, the percentage of EBs that corre- 
late with downflows in area 1 is 86%, as mentioned 
above. Local minima of upward Doppler velocity 
are seen for 6% of EBs in area 2 and for 1 % of 
EBs in area 3. The majority of EBs in these areas, 
however, are associated with actual downflows. 

An interesting feature of EBs in area 3 is that 
they occur in unipolar areas. These EBs are trig- 
gered above the leading edges of rapidly mowing, 
recently emerged magnetic flux bundles. We con- 
clude that these EBs occur in locations of int erec- 
tion between the pre-existing magnetic field lines 
and the newly emerged magnetic flux. As we will 
show in § 4.4, the interacting magnetic config- 
urations have different geometries. Parts oif the 
pre-existing magnetic field are relatively horizon- 
tal, while the moving magnetic patches form large 
angles with the horizontal plane. The proper - mo- 
tions of these patches force the different magnetic 
configurations to interact, producing swarms of 
EBs. Examples are given in § 4.4 (Figs. S and 
7) and we attempt to interpret this interesting re- 
sult in § 7. 

4.3. EB association with 1600 A BPs in the 
temperature-minimum region 

In Severny (1968), and in references therein, 
EBs are classified as photospheric entities. Re- 
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Magnetic neu- Proper motions 
tral fines of EB loci 

Downflows 

1600 A BPs 

Area 1 

... 77% 

(< 0.2 — 0.4) km s ~ 1 

86% 

54% 



(supergranular boundary) 

(—0.4, —0.1) km 




(0.4 — > 0.8) km s~ x 





(cell interior/MDFs) 



Area 2 

... 65% 

(< 0.2 - 0.4) km s" 1 

86% 

56% 




(—0.3, —0.1) km 

s -1 

Area 3 

... 36% 

(< 0.2 — 0.8) km s _1 

81% 

55% 




~ —0.5 km s~ l 



Table 1: Synoptic results on the association of EBs with the underlying photospheric pattern (magnetic fields 
and mass flows) and the association of EBs with 1600 A BPs in the T min region. 


cently, however, it has been conjectured that EB 
emission originates from close to the T min , at a 
height 600 km - 1100 km above the 75000 = 1 
layer (Payne 1993). BPs in 1600 A are a typi- 
cal feature of the T min , occurring due to the ion- 
ization of Si and Fe in wavelengths 1500 - 1680 
A (Vemazza, Avrett, & Loser 1981). Qiu et al. 
(2000) reported that ~ 61% of their EBs correlate 
with 1600 A BPs. 

Rrom the 13 TRACE frames at 1600 A that are 
nearly simultaneous with FGE Ha filtergrams, we 
■visually examined 448 EBs (142 EBs in area 1, 106 
EBs in area 2, and 200 EBs in area 3) and their 
association with 1600 A BPs. In agreement with 
Qiu et al. (2000), we found that slightly more 
than half EBs are associated with 1600 A BPs, in- 
dependently of the area. In 54%, 56%, and 55% of 
EBs from area 1, area 2, and area 3, respectively, a 
clear correlation with 1600 A BPs could be seen. 
These percentages should be considered average, 
however, since the association of EBs with BPs 
depends on the size of the EBs. For bright and 
persistent EBs the percentage exceeds 90%, but it 
can be as low as 45% for faint, short-lived EBs. 
We noticed an interesting property of EBs asso- 
ciated with MDFs in the supergranular interior. 
Although these EBs typically persist for as long 
as the MDFs me visible, no 1600 A BPs appear 
related to them. 

Our results are in agreement with the results 
of Qiu et al. (2000). An interesting open prob- 
lem is to determine which emission (Ha or UV) 
precedes the other. Unfortunately, our moderate 
cadence does not allow such an investigation. The 


photospheric conditions beneath EBs, as well as 
the association of EBs with 1600 A BPs, are sum- 
marized in Table 1. 

4.4. Some representative EB cases 

In Fig. 5 we depict three EB examples (one for 
each of the areas 1, 2, and 3). Ha — 0.8A emission 
is outlined by the thick, solid contours, "while 
the thick dashed contours outline 1600 A BPs. 
The grayscale background corresponds to both the 
LOS longitudinal magnetic field (left column; im- 
ages labeled (a)) and the Doppler velocity (right 
column; images labeled (b)). Dark patches corre- 
spond to vectors oriented downwards, while bright 
patches correspond to vectors pointing upwards. 
The thin solid contours represent neutral magnetic 
lines (left column) and zero-velocity lines (right 
column). 

In all cases, EBs are clearly associated "with 
1600 A BPs. The EB from area 1 (upper row) oc- 
curs above a well-defined neutral line and an. area 
of photospheric downflows. The EB from area 2 
(middle row) also correlates with a magnetic neu- 
tral line. This line is more noisy compared to 
the one shown in the upper row, since magnetic 
fields in area 2 are typically weaker compared to 
those found in area 1. The EB is, again, associated 
with photospheric downflows. The EB from, area 
3 (lower row) clearly occurs in the absence of neu- 
tral magnetic lines, in a unipolar area of positive 
polarity. Nevertheless, this EB is also associated 
with photospheric downflows. 

In Fig. 6 we show two examples of EBs for each 
of the areas 1, 2, and 3. We have now defined 
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Fig. 5.— Three examples of EBs and their local conditions. The locations of enhanced Ha - 0.8 A pmissinn 
axe shown by thick solid contours. Enhanced UV 1600 A continuum emission is indicated by dashed contours. 
The grayscale backgrounds correspond to the (nearly) simultaneous LOS longitudinal magnetic field (left 
column; cases labeled (a)) and Doppler velocity (right column; cases labeled (b)). Each EB is plotted twice 
m each row using the above-mentioned backgrounds. The thin solid line indicates neutral magnetic lines 
(left column) and zero-velocity lines (right column). The tick mark separation is 1 arcsec. The UT times 
of the observations axe shown at the lower left comer of the images, for the Ha — 0.8 A filtergrams (FGE) 
and the UV images (TRACE). The exact locations of the EBs in the AR are indicated by arrows in Fig. 2. 


a number of one-dimensional spatial profiles 
(line-cuts) and we depict the evolution of various 
parameters along the selected line segments. In- 
formation includes the Ha - 0.8 A contrast (thick 
solid fine), the 1600 A contrast (dash-dotted line), 
the LOS longitudinal magnetic field (thin solid 
line), and the Doppler velocity (dashed line). The 
background for the Ha and the UV emission i s set 
to zero. 

The EB in Fig. 6a is associated with both a 
well-defined neutral magnetic line (separating 


magnetic fields ~ 500 G and ~ -1000 G) and 
a 1600 A BP. The Doppler motions are ascend- 
ing. However, there is a local velocity minimum 
of ~ 0.1 km s -1 roughly coinciding with the EB, 
whereas the Doppler velocity in the surround- 
ing areas is as high as 1 fan r 1 . In Fig. 6b, 
the brighter EB is associated with a downflow 
~ —0.5 km s ~ 1 , a neutral magnetic line and a 
1600 A BP. In area 2 (Figs. 6c, d), both the mag- 
netic fields and the Doppler motions are weaker, 
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Fig. 6. — Six examples of one-dimensional spatial profiles of EBs. The information along the line is the 
Ha- 0.8 A contrast (thick solid curves), the UV 1600 A contrast (dash-dotted curves), the LOS longitudinal 
magnetic field (thin solid curves) and the Doppler velocity (dashed curves). A legend is given in diagram 
(d), for reference. Each curve is the average of 10 similar curves obtained from respective line-cuts. The 
line-cuts are parallel and uniformly distributed. The distance between adjacent line-cuts is equal to the pixel 
size. 


compared to those of area 1. In Fig. 6c, a rather 
faint neutral line exists beneath the EB. The 
Doppler velocities, however, show a more profound 
behavior, reaching values ~ —0.15 km s~ x be- 
neath the EB. At the EB’s surroundings, Doppler 
velocities are of the order 0.4 km s" 1 . No neu- 
tral line can be seen in Fig. 6d, for both EBs. 
The fainter EB is not associated with a down- 
flow or a 1600 A BP, but a strong downflow of 
~ —0.6 km 1 and a 1600 A BP correlate nicely 
with the brighter EB. In Figs. 6e, f (area 3), the 
longitudinal magnetic field is clearly unipolar. A 
downflow of the order —0.1 km s -1 , however, 


occurs beneath the EB of Fig. 6e, while the sur- 
roundings of the EB show upflows of the order 
(0.12,0.24) km s -1 . The EB coincides precisely 
with a UV BP. In Fig. 6f one notices two 1600 A 
BPs but only one EB. The EB occurs above the 
edge of an area dominated by descending motions 
of the order —0.5 km s '" 1 . 

In Figs. 7a, b we use the same line-cuts as in 
Figs. 6e, f, respectively. The information shown 
here, however, is the Ha — 0.8 A contrast (solid 
lines) and the zenith angle 0 of the photospheric 
magnetic field vector with respect to the LOS di- 


15 









o 2 4 6 8 10 12 

Space (arcsec) 



Space (arcsec) 

Fig- 7. — One-dimensional spatial profiles for the 
EBs of the last two examples of Fig. 6. The EBs 
from Figs. 6e and 6f are the same with the EBs of 
the diagrams (a) and (b), respectively. The infor- 
mation along the line is the Ha - 0.8 A contrast 
(solid curves) and the zenith angle of the mag- 
netic field vector (dashed curves). The orienta- 
tion of motion of the relatively vertical magnetic 
flux bundles against the relatively horizontal, pre- 
existing field, is from right to left. 

rection (dashed lines). 6 ranges between 0° (the 
magnetic field vector points outwards, along the 
LOS axis) and 180° (the magnetic field vector 
points inwards). For 0 ~ 90° the magnetic field 
vector is nearly transverse to the LOS axis. From 
Fig. 7 we find that there is a steep gradient of the 
zenith angle in the locations of EBs. The newly 
emerged magnetic patches which drift toward the 
leading sunspot in area 3 are mostly perpendicu- 
lar to the photosphere (25° < 0 < 30°) and they 
interact with the mostly horizontal, pre-existing 
magnetic fields (70° < 9 < 80°). Therefore, in 
both examples of Fig. 7, EBs occur where rela- 
tively vertical magnetic fields are interacting with 
pre-existing, relatively horizontal, fields. 


5. The statistical properties of EBs 
5.1. The fractal character of EB clusters 

Let us now investigate statistically the spatial 
distribution of EBs. A first-order approach to 
tackling this problem is the investigation of the 
fractal character of EB occurrence. For large-scale 
structures, the numerical calculation of the frac- 
tal dimension provides information on the degree 
of complexity that characterizes the structures’ 
boundaries. For small-scale structures, occurring 
in swarms and having sizes comparable to the spa- 
tial resolution, the fractal dimension provides in- 
sight on how inhomogeneous the clustering pro- 
cess of these structures is. This information, will 
be provided for EBs, since EBs meet the second 
type of conditions. 

To calculate the fractal dimension, we rely on 
a two-dimensional box-counting method (for a 
complete description see, e.g., Mandelbrot 1983 
and references therein) as follows: We cover each 
Ho. — 0.8 A frame with a two-dimensional, uni- 
form, rectangular grid, consisting of boxes with 
size A x A. If our field of view has a size L x L, 
then it is covered by [LJ A) 2 boxes with a non- 
dimensional area £ 2 ; e = A. The total number of 
boxes being [LJ A) 2 , we count the number of boxes 
that enclose non-zero information, that is, at least 
one EB or a part of at least one EB’s boundaries. 
If this number is N( A), we vary the box-size A (or, 
equivalently, the non-dimensional box-size e) and 
study the variability of N( A) vs. A (or A). For a 
fractal structure, one should obtain 

N(e) oc (i) c » (1) 

where T> 0 is the (box-counting) fractal dimension. 
T>o can be at most equal to the embedded (Eu- 
clidean) dimension of the space, for non-fi-actal 
structures, or smaller, for fractal structures , but 
not larger than that. In our case, V 0 < 2. P 0 = 2 
would imply a spatially homogeneous event distri- 
bution, whereas for T> 0 < 2 (and large statistical 
samples) events show a strong preference to orga- 
nize in clusters. 

An example of the numerical calculation of the 
fractal dimension is shown in Fig. 8a. Moreover, 
we have calculated the fractal dimen sio n of EB 
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Background contrast 


Fig. 8. — (a) Numerical calculation of the fractal 
dimension of EB clusters for a given frame (shown 
in the inset). The dynamical range over which the 
fractal dimension is estimated exceeds one order 
of magnitude, (b) The time variability of the frac- 
tal dimension of EBs for I Co = 0.08. The mean 
value of the fractal dimension (~ 1.4) is shown 
by the dashed line, (c) The mean fractal dimen- 
sions for different background selections. The er- 
ror bars correspond to the standard deviation for 
each fractal dimension, as calculated by the square 
root of the variance in the time variability of the 
fractal dimension for a given background. Clearly, 
the standard deviation (or the variance) increases 
with increasing background. For I Co > 0.11 the 
calculation of the mean fractal dimension becomes 
unreliable, due to poor statistics. 


clusters for all FGE Ha - 0.8 A filtergrams, thus 
constructing the timeseries of the fractal dimen- 
sion variations for different background contrasts 
in the range (0.05,0.20). We found that EB clus- 
tering is a fractal process, since we obtain a fractal 
dimension Vo between 1.5 (for I Co = 0.05) and 1.0 
(for I Co = 0.20). Moreover, the fractal dimension 
does not change appreciably with time for a given 


background, provided that the number of EBs is 
large enough to allow a reliable calculation. Stated 
differently, the fractal dimension fluctuates around 
a mean value. We show this behavior in Fig. 8b, 
for I CQ = 0.08. The mean fractal dimension is 
V 0 — 1.4. Fig. 8c shows the average (over the 
observing time) calculated fractal dimens ions for 
different backgrounds I Co . We find that the fractal 
dimension decreases with increasing 7 Co . 

The above findings suggest, from a different 
perspective, that the spatial distribution of EBs is 
inhomogeneous and that preferential locations in 
the AR are prolific in EB production. As we know, 
these locations are the supergranular boundaries 
(area 1) and area 3, where new magnetic flux ac- 
cumulates on the leading sunspot (Fig. 2). Sec- 
ondary locations are the polarity reversal area and 
the supergranular interior. The relative constancy 
of the fractal dimension with time (Fig. 8a) in- 
dicates that the character of the clustering pro- 
cess does not change significantly within our ob- 
serving interval. Changes in the fractal dimension 
with time, for a given background contrast, can 
be mostly attributed to the varying numbers of 
EBs in each filtergram. This behavior is, proba- 
bly, expected because (i) the preferential locations 
of EB occurrence do not change within our observ- 
ing time, and (ii) the timescale of flux emergence, 
which ultimately triggers EBs, is large compared 
to our observing interval. The decreasing mpan 
fractal dimension with increasing background con- 
trast is due to the significantly smaller EB num- 
bers as the background increases. For as long the 
statistics are sufficient, i.e., for I CQ < 0.11 (Fig. 
8b), this result indicates an enhancement of the 
clustering process. For I Co >0.11, the rapidly de- 
creasing numbers of EBs make the calculation of 
the fractal dimension unreliable, as shown t>y the 
large error bars in Fig. 8b. 

5.2. Definitions and statistical properties 
of typical EB parameters 

In this analysis, we define an array of parame- 
ters that highlight different aspects of EB “size” . 
The parameters to be defined are the area of EBs, 
the EB lifetime (duration), the maximum and the 
total integrated EB contrasts. 

To account for finite resolution effects, since 
EBs have sizes close to the spatial resolution, we 
introduce the following correction: We assume 
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that the peak contrast I p (t) of an EB, at a given 
frame (time) t, decreases as a gaussian function 
of the distance l from the location of the peak. 
The parameters of the gaussian may be then cal- 
culated, using the EB length /(f), as deduced by 
the measured EB area, for a given background con- 
trast I Co (Appendix). Then we calculate the cor- 
rected area 5(f) for which the intensity fades to 
the mean intensity J 0 (contrast equal to 0). This 
area is, thereafter, treated as the EB area at time 
t. 

This being said, we define the following EB pa- 
rameters: 

(1) The lifetime D of EBs. This is estimated 
simply by the number of frames over which 
an EB persists, multiplied by the average 
elapsed time between frames (~ 7.5 min). 
We mentioned before that a calculated EB 
lifetime may be the cumulative result of 
more than one EBs that recur and that 
short-lived EBs may have been unaccounted 
for. For this reason, we will not pursue any 
robust conclusion based on EB lifetimes. We 
have chosen this parameter, however, be- 
cause it is required for the definitions of the 
parameters that follow. The duration of EBs 
is the only EB parameter that remains un- 
affected by the gaussian correction. 

(2) The corrected maximum EB area S max . This 
parameter is provided by the mrimnm 


where I(S,t) is the contrast of the EB emis- 
sion throughout the area 5(f). Applying the 
gaussian correction, we obtain the corrected 
integrated contrast 

I s (t) ~ I p {t)b^ erf (0.745 1|/2 (t) /26) (5) 

where 6 2 = erf(x) = e~^du 

is the error function of x, I p (t) is the mea- 
sured peak EB contrast, and S(t) is the cor- 
rected EB area, at time t. For the derivation 
of equation (5), see Appendix. 

The corrected maximum integrated contrast 
is now given by the maximum 

Imax = max{I(t '); t' = t, t + 1, ..., 1 +• D} 

( 6 ) 

of the corrected integrated contrasts over the 
lifetime D of an EB. The maximum inte- 
grated contrast of an EB refers to the EB’s 
peak activity and is insensitive to the recur- 
rence effect. If multiple, co-spatial EBs, are 
detected as one event due to undetected re- 
currence, then I max reflects the peak activity 
of the most prominent of them. 

(4) The corrected total integrated contrast I tot of 
EBs. This parameter occurs after integrat- 
ing in time the spatially integrated contrast 
I s (t) of EBs (eq. [5]) for the EBs’ lifetime 
D, namely, 


Smax = max{S(t '); t' = t,t + 1, ..., t + D] 

( 2 ) 

of the corrected areas S(t‘) over an EB’s life- 
time D . The corrected area S(t) of an EB 
at time (frame) t is given by 


S(t) = S p (t) 


l°9lp{t) 

l°9{Ip{t)/Ic o) 


( 3 ) 


with Ip(i) and S p (t) being the maximum EB 
contrast and the measured EB area at time 
t. The derivation of equation (3) is provided 
in detail in the Appendix. 


(3) The corrected maximum integrated contrast 
Imax of EBs. The spatially integrated con- 
trast I s (t) of an EB, at a given time £, is 
provided by the integral 

*•(*)=/ I(S,t)dS (4) 
Js(t) 


rt’+D 

Itot = J (7) 

The total integrated contrast is a measure of 
the total activity of EBs and is, obviously, 
affected by the recurrence effect. 

The analysis hereafter is totally based on param- 
eters obtained under the gaussian correction for 
unresolved effects. 

In Figs. 9a, 9c, 9e, and 9g we provide th.e dif- 
ferential distribution functions of the duration, the 
maximum area, the maximum integrated contrast 
and the total integrated contrast of EBs, respec- 
tively, using I CQ = 0.05 as the background con- 
trast. 

The central feature of the distribution functions 
of Fig. 9 is their power-law shape. Power laws are 
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Fig. 9. Distribution functions and power-law index variations for typical EB parameters: Duration dis- 
tribution (a) and duration-index variation (b), maximum area distribution (c) and index variation (d), 
maximum integrated contrast distribution (e) and index variation (f), and total integrated contrast distri- 
bution (g) and index variation (h). The best least-squares power-law fits are shown by straight lines (left 
column). Error bars (right column) show the standard deviation calculated for each index. The mean value 
of the calculated power-law indices is shown by dashed straight lines (right column). 


suggestive of the existence of an intrinsic self- 
similarity in the triggering process of EBs. Self- 
similarity, or scale-invariance, in the distribution 
functions of EB parameters fits well with the frac- 
tal properties of EB clustering. This core property 
of EBs suggests that the same behavior should be 
expected for undetected EBs, as well, down to the 
elementary EB sizes. Power laws in Fig. 9 are of 


the form 


dN(X) 

dX 


<xX~ Qx 


( 8 ) 


where X — {-D, S max , I max > ^tot} and cnx is the 
respective power-law index. We obtain power laws 
for all typical EB parameters, independently of the 
background selection, provided that the EB 
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Fig. 10. Correlation scatter plots linking typical EB parameters for 7 Co = 0.05. The straight lines 
correspond to the least-squares power-law best fit. The scatter plots correlate the maximum integrated 
contrast and the total integrated contrast (a), the maximum integrated contrast and the maximum area 
(b), the total integrated contrast and the maximum area (c), the maximum area and the duration (d) the 
maximum integrated contrast and the duration (e), and the total integrated contrast and the duration (f). 
The straight-line components seen in plots (a) and (c) correspond to EBs with lifetime equal to 1 (frame), for 
which I tot — I max- These events have been excluded from the fitting process. The quality of the correlations 
is measured by means of the linear correlation coeflicient (cc). 


numbers are sufficient. In Figs. 9b, 9d, 9f, and 
9h we show the variations of the power-law indices 
ax for the distribution functions of Figs. 9a, 9c, 9e 
and 9g, respectively, for different backgrounds / Co . 
The error bars in the calculation of the power-law 
indices become larger as the background increases, 
since the EB sample becomes smaller and smaller. 

We stated previously that we will not pursue 
conclusions based on the EB duration, because of 
our moderate cadence. Indeed, from Fig. 9a 


one notices durations well above 100 min, which 
clearly suggest that the duration distribution func- 
tion is severely affected by the recurrence effect. 
Similar EB durations have been obtained in early, 
low-cadence studies (Bruzek 1972; Vorpahl and 
Pope 1972). 

Although power laws themselves should be 
viewed as the main result of Fig. 9, it is note- 
worthy to discuss briefly the values of the power- 
law indices. We made tests without applying the 
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gaussian correction and we obtained significantly 
flatter power laws compared to their respective 
ones in Fig. 9. The correction, therefore, did not 
change the power-law shape of the distribution 
functions, but it resulted in different power-law in- 
dices. It should be expected that different adopted 
assumptions give different power laws in all distri- 
bution functions, except the duration distribution 
function which is not subject to any correction. 
The mean power-law index ab for EB durations 
is ab — 2.20. The mean power-law indices for the 
maximum integrated contrast and the total inte- 
grated contrast are aj m ~ 3.53 and aj t ~ 2.07, 
respectively. The mean power-law index a s of the 
EB areas is as — 2.44. 

We have, furthermore, investigated possible 
correlations between typical EB parameters. To 
achieve this, we have constructed correlation scat- 
ter plots linking these parameters. These scatter 
plots are shown in Fig. 10. Again, we obtain 
power-law relations between EB parameters, of 
the form 

X oc Y 5xy (9) 

where X,Y == {D, S max ,I max , Itot} and Sxy is 
the correlation index between parameters X and 
y (&xy = 1 for X = cY; c constant). The (lin- 
ear) correlation coefficients are moderate-to-high, 
varying in the range (0.55,0.99). Notice the ex- 
cellent relation between maximum EB areas and 
maximum integrated contrasts (Fig. 10b). A high 
correlation coefficient should be expected, since 
the maximum integrated contrast of an EB, in 
most cases, occurs when the EB area peaks. This 
striking correlation, however, points to an invari- 
ant relation between the area and the integrated 
contrast of EBs, independently of the background 
selection. This result is not due to the gaussian 
correction; it is also obtained when the correction 
is not applied. 

The above is the first statistical study of sev- 
eral EB parameters. The power-law distribution 
functions and the interrelated event parameters, 
however, are well documented results in the sta- 
tistical studies of flares, soft X-ray, and extreme 
ultraviolet (EUV) BPs. Remarkably, the quality 
of the correlations is even better for EBs, in some 
cases. We will discuss this result in § 8. 


6. Energetics of EBs 

We will now estimate typical EB energies. En- 
ergy losses in the EB volume, as well as mo- 
tions associated with EBs, will serve as proxies 
of the actual released energies. Our results w r ill be 
compared to previous estimations of EB energies, 
which vary in the range (10 25 ,10 27 ) erg (Bruzek 
1972; Hu, Song, & Li 1995; Henoux, Fang, & Ding 
1998). 

Energy losses during EBs should be attributed 
to radiation and thermal conduction. The con- 
ductive losses may be neglected in the low chro- 
mosphere, because of the small temperature gra- 
dients (Payne 1993). Therefore, the dominant 
energy loss mechanism in EBs is radiative cool- 
ing. We adopt the expression of Nagai (1980) for 
the net radiative loss rate per unit volume e ra d 
(erg cm" 3 s'" 1 ), namely, 

e r ad ^ a(T)n 2 xg(T) (10) 

where a(T) is the radiative reduction coefficient, 
n = n e + uh — n e , with n e , n# the numeri- 
cal densities of electrons and neutral hydrogen, 
respectively, x is the ionization degree and g(T) 
is a semi-empirical function of the temperature 
T. Adopting Nagai’s (1980) findings, we obtain 
g(T) ~ 10 ~ 23 erg cm 3 s _1 , for T = 10 4 AT, as- 
suming that this temperature is typical of the EB 
volume. For T < 5 x 10 4 K, moreover, the ion- 
ization degree is approximated by the expression 
(Brown 1973) 

, X 2 , 51405 

1o 9yZ ^ = ~ l °9 n f~ + 0.5 logT + 16.0039 

( 11 ) 

Payne (1993) used a numerical density n ~ n e — 
2.24 x 10 12 cm -3 in the EB volume. This den- 
sity is consistent with works on radiative txans- 
fer models (Canfield and Gayley 1987; Hawley 
and Fisher 1994; Abbett and Hawley 1999) of 
solar flare atmospheres. These models predict a 
broadening of the Ha line, similar to that ob- 
tained in EBs, when a sharp temperature increase 
(AT ~ 10 3 - 10 4 K) occurs, followed by an abrupt 
increase of the ionization degree. The high ioniza- 
tion degree allows plasma confinement due to the 
fact that the magnetic field is stronger than the 
equipartition value and the density is kept in the 
range (10 12 ,10 13 ) cm -3 . Adopting n ~ 10 12 cm~ 3 
and T — 10 4 K in EBs, we obtain, from equation 
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(11). X — 0.89. Moreover, we obtain (Nagai 1980) 
a(T) — 0.1 for T ~ 10 4 K. Therefore, we find 
from equation (10) that the net radiative loss rate 
in EBs is £ r ad — 0.89 erg cm~ 3 s -1 . 

Knowledge of the EB radiative volume and the 
EB duration are required to calculate the peak 
and the total radiative losses during EBs. The EB 
volume V E b may be approximated by the relation 
Veb ~ Sh, where S is the area occupied by an 
EB and h is the vertical extent (height) of EBs. 
We have established the association of EBs with 
1600 A BPs in the T m ; n region which allows a 
conservative estimate of the EB height, that is, 
h ~ 10 2 km. Other authors (Kitai and Muller 
1984) argue that even a height h ~ 10 3 km should 
not be ruled out. 

The peak radiative loss rate (peak luminosity) 
Prod (erg s -1 ) of EBs is given by the relation 

P rod = £radfVEB(max) (12) 

where V EB(max) ~ S max h is the maximum volume 
attained by a given EB. This volume is assumed 
linearly related to the maximum area S max , calcu- 
lated for this EB, and a typical height h = 10 2 km. 
f is the radiative filling factor, which is consid- 
ered unity, thus assuming that the entire EB vol- 
ume radiates during the energy release process. 
To calculate the total radiative losses E rad over 
the EBs’ lifetime, we further need to adopt a typ- 
ical temporal profile of EBs: EBs are assumed ei- 
ther symmetric with a plateau roughly equal to the 
peak released energy (Payne 1993) or symmetric 
with roughly equal rise and decay times (Qiu et 
al. 2000). We adopt the second assumption, thus 
obtaining 

Er ad ~ (13) 

where D is the lifetime of an EB showing a peak 
luminosity equal to Prod- The first assumption 
provides a larger E rad , by a factor of two. Al- 
though we have roughly estimated the durations 
of EBs in our sample, we underlined that our es- 
timations should not be considered reliable, due 
to the moderate cadence, compared to typical EB 
lifetimes, and the recurrence effect. Indeed, the 
latter effect is obvious in Fig. 9a. Therefore, we 
will use typical EB lifetimes, provided by previous 
authors. The adopted typical value of EB lifetimes 
is Deb — 600 s (10 min). We have, moreover, ap- 



Peok radiative loss rate P™, (erg s' 1 ) 

Fig. 11. — Differential distribution function of the 
peak and the total radiative losses in EBs for a 
background contrast I Co = 0.05. Calculations have 
been performed assuming a vertical extent of the 
EB volume equal to 100 km, an EB duration equal 
to 10 min and a temperature of the radiative vol- 
ume equal to 10 4 K. The straight line is the least- 
squares power-law best fit. 

plied the correction for unresolved effects to the 
calculation of the maximum EB areas S max . 

The differential distribution functions of the 
peak and the total radiative losses in EBs are given 
in Fig. 11 for I Co = 0.05. The peak radiative losses 
in EBs vary in the range (4 x 10 24 , 10 26 ) erg s -1 , 
while the total energy losses due to radiative cool- 
ing of EBs vary between ~ 1.2 x 10 27 erg and 
~3x 10 28 er<7. Another notable feature in Fig. 11 
is the power-law shape of the distribution function. 
The (peak and total) radiative energy losses E of 
EBs obey a distribution function oc JE~ a , 

with power-law index a = 2.12 ± 0.05. It should 
be noted that the power law obtained in Fig. 11 
depends largely on the power law of the maximum 
EB areas (Fig. 9c) and, therefore, on the gaussian 
correction. Moreover, it depends on the fact that 
we adopt a constant, typical EB lifetime. Conclu- 
sions relying on the value of the power-law index 
should be considered with care. 

On the other hand, EBs are often associated 
with kinematic events. Rust (1968) reported the 
correlation of EBs with surges, ejecting witli ve- 
locities ~ 10 2 km s ~ 1 and showing threads ~ 
5 x 10 3 km high. Assuming a thread-like surge 
with a height ~ 5 x 10 8 cm, an aspect ratio of, say, 
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10, that is, a diameter of ~ 5 x 10 7 cm (500 km ), 
an electron density ~ 10 n cm~~ 3 and a velocity 
~ 10 7 cm s ~ x , then the kinetic energy of the surge 
is ~ 4.5 x 10 26 er#. The radiative energy of the 
surge is ~ 10 26 er#, thus giving a total surge en- 
ergy ~ 5.5 x 10 26 er#. Although the association of 
EBs with surges is not clear in this work, we find 
that a total surge energy of ~ 5.5 x 10 26 er# is com- 
parable to the total radiative losses in EBs, at least 
in case of the weaker events. If the released energy 
during EBs triggers surges then, based on our cal- 
culations, the total radiative losses and kinematic 
effects axe balanced by a total released energy in 
the range (1.8 x 10 27 ,3 x 10 28 ) erg. This estima- 
tion is, in general, consistent with the previously 
reported energy range of EBs, although it is some- 
what shifted toward higher energies. 

If the numerical density in the EB volume is of 
the order 10 12 cm -3 , what is the temperature en- 
hancement in EBs ? For the net radiative energy 
loss rate e ra d to be balanced by a thermal energy 
input, we require a temperature enhancement AT, 
given by the relation 

£rad - \nk B AT or AT ~ (14) 

2 Ztikb 

where n ~ n e + rii ~ 2 x 10 12 cm -3 and 
ks is the Boltzmann constant. For € rad = 
0.89 erg cm“ 3 s _1 , we obtain, from equation (14), 
AT ~ 2.15 x 10 3 K. This temperature enhance- 
ment is not beyond reason and, in fact, it is close 
to the temperature enhancement of ~ 1.5 x 10 3 K, 
suggested by Kitai (1983). To account for both 
radiative losses and mass motions in surges, a well- 
defined EB with S m ax ~ 10 arcse(? and duration 
D ~ 10 min requires a temperature enhancement 
AT ~ 4.3 x 10 3 K. In the above analysis, we 
have implicitly neglected the possibility of a non- 
thermal energy release component which may, as 
well, be present during EBs. Based on our thermal 
calculations, the rate of energy release, for a rea- 
sonable AT ~ 2.15 x 10 3 AT, is comparable to the 
radiative cooling rate. The cooling timescale re- 
sponding to a single energy deposition, therefore, 
should be short, of the order of a few s . This result 
is consistent with the results of non-LTE radiative 
transfer models (Abbett and Hawley 1999) that 
predict short cooling timescales in heated plas- 
mas in the chromosphere, characterized by a line 
broadening in Ha. 


7. An interpretation of EB triggering 

Our results (§ 4) suggest that EBs occur on the 
interface of colliding magnetic configurations, ei- 
ther above or in the absence of photospheric mag- 
netic neutral lines. An intriguing feature of such 
complex topologies is the likely presence of qiiasi- 
separatrix layers (QSLs - van Driel Gesztelyi et 
al. 1994; Demoulin et ai. 1996; Demoulin et aJ. 
1997). QSLs are a generalization of separatrices, 
which apply to magnetic null points, to tliree- 
dimensional configurations which do not neces- 
sarily require nulls (Demoulin and Priest 1997). 
QSLs, or separatrices, are formed on the interface 
between different flux systems, that is, where the 
magnetic connectivity pattern changes. It is ex- 
pected, therefore, that separatrices and QSLs are 
formed (i) on the dipolar supergranular bound- 
aries, (ii) on MDFs, and (Hi) between emerging 
vertical fields and pre-existing horizontal fields at 
the vicinity of the leading spot, since the two 
topologies are geometrically different. The stbove 
areas are all prolific in EB production. It is 
likely, therefore, that EBs occur on separatrices 
and QSLs in the low chromosphere. 

In Fig. 12 we show schematically how our data 
support this triggering cause for EBs. In Fig. 12a 
we show a magnetic configuration similar to the 
one expected for the supergranular boundaries. 
The strong neutral lines are sustained by converg- 
ing flows and an X-point separatrix forms between 
the two flux systems. The EB occurs on the sep- 
aratrix, above the photospheric neutral line. In 
Fig. 12b, small X-point separatrices are formed 
above MDFs. The MDF configurations in both 
our data and Fig. 12b appear inclined toward the 
direction of the flow. The EBs occur above the 
neutral lines of the MDFs. Our MDFs exhibit a 
U-loop shape, instead of the classical, emerging, 
fl-loops (Bernasconi et al. 2002). The absence of 
apparent links (magnetic field lines) with higher 
atmospheric layers probably explains the latck of 
correlations between EBs above MDFs and 1600 
A BPs, pointed out in § 4.3. In Fig. 12c, a QSL is 
formed between the pre-existing, nearly horizon- 
tal magnetic field and the relatively perpendicular 
moving magnetic flux bundles. In our data, this is 
shown by the gradient of the zenith angle between 
the two configurations in area 3 (Fig. 7). The EB 
of Fig. 12c occurs in the unipolar area of the inter- 
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Fig. 12. — Three different mechanisms of EB trig- 
gering on separatrices and QSLs. (a) Magnetic re- 
connection in strong, dipolar magnetic fields sus- 
tained by converging horizontal flows, (b) EBs 
triggered by magnetic reconnection above MDFs. 
(c) EB triggered by the interaction of two topolog- 
ically different, unipolar, magnetic configurations. 


acting, topologically different magnetic structures. 
Our results regarding the magnetic field topology 
in area 3 are consistent with the results of bites, 
Skumanich, & Martinez Pillet (1998). The rela- 
tively horizontal fields form the low-lying magnetic 
canopy with weak magnetic fields, that connects 
the two sunspots. The relatively vertical fields are 
the newly emerged individual flux elements that 
attain significant horizontal velocities and strong, 
vertical magnetic fields when they migrate away 
from the flux emergence zone. In all shown sce- 
narios of EB triggering, in Fig. 12, EBs occur 
close to the photospheric footpoints of magnetic 
field lines. This probably explains the downflows 
associated with almost all EBs, as being drainage 
motions of previously elevated, cold material, to- 
ward the photospheric footpoints of emerged flux 
ropes (Bruzek 1967; Zwaan 1985). 

Separatrices and QSLs (hereafter QSLs, for 
brevity) are candidate areas for magnetic recon- 
nection on the photosphere/low chromosphere. In 
fact, they may be the only locations in the lower 
atmosphere that allow a large growth rate of the 
tearing mode instability, especially at the vicinity 
of the T m i n region (Sturrock 1999). Magnetic 
reconnection normally requires current sheets 
formed by anti-parallel magnetic fields but we ob- 
tain EBs in both dipolar and unipolar areas. How- 
ever, in both theory (Soward 1982) and numeri- 
cal simulations (Linton, Dahlburg, & Antiochos 
2001), it has been reported that reconnection can 
be efficient even at small contact angles between 
the interacting magnetic field lines. Free mag- 
netic energy is accumulated on QSLs. Stochastic 
magnetic reconnection in a turbulent environment 
effectively releases the free energy and fragments 
QSLs. Our data favor EB triggering by magnetic 
reconnection on QSLs, due to the turbulent evolu- 
tion of the magnetic field. This process probably 
accounts for (i) the self-similarity in EB statistics, 
(ii) the small sizes of EBs, and (Hi) the preferred 
elongated EB shape, reminiscent of a Petschek 
reconnection scenario. 

8. Summary and Discussion 

EBs occur in areas of strong, colliding magnetic 
configurations. In particular, they occur: 

(1) Above magnetic neutral lines outlining the 
boundaries of supergranular cells. 
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(2) Above the magnetic neutral lines of MDFs. 

(3) In the absence of neutral lines, when topo- 
logically different magnetic configurations 
interact. In this case, EBs are triggered due 
to the interaction between the pre-existing, 
low-lying horizontal magnetic canopy and 
plage-like, strong vertical fields, at the vicin- 
ity of the leading sunspot. 

The above locations are most prominent in EB 
triggering throughout our entire observing time, 
thus resulting in an inhomogeneous spatial distri- 
bution of EBs. EBs appear in clusters, or swarms, 
in these preferential locations. The recurrence ef- 
fect of EBs originates from this behavior. 

The correlation of EBs with photospheric mass 
flows is as follows: 

(1) The Doppler photospheric velocities beneath 
EBs indicate downflows. In several cases, 
EBs occur above the edge of areas domi- 
nated by downflows. We have seen EBs asso- 
ciated with upflows but, in these cases, there 
is a relative downflow in the locations of 
EBs, with respect to the EBs’ surroundings. 
Downflows probably occur due to drainage 
motions of cold material toward the photo- 
spheric footpoints of previously emerged flux 
ropes. 

(2) EB loci move parallel to transverse photo- 
spheric flows, as calculated by correlation- 
tracking techniques applied on the FGE 
white-light images. This behavior is clearly 
seen in areas of systematic transverse flows, 
such as the supergranular area and the vicin- 
ity of the leading sunspot. On the bound- 
aries of the supergranule, which are sus- 
tained by converging transverse flows, EB 
loci move with smaller velocities. 

The typical EB size is 1.8 x 1.1 arcsec?. EBs 
exhibit an elongated shape with their length 40% 
larger than their width. A large number of EBs 
have sizes close to, or even below, our spatial res- 
olution (0.8 arcsec). This suggests that improved 
spatial resolution may set the typical EB sizes to 
even lower values. We found a mean EB occur- 
rence rate ~ 1.43 EBs min” 1 for an EB-producing 
area ~ 1800 arcsec ? . Given our moderate cadence 


(7.5 min), the recurrence effect, our finite spa- 
tial resolution, and the overlaying AFS filaments 
which may obscure several EBs, we conclude that 
this occurrence rate is only a lower limit of the 
actual occurrence rate of EBs. 

The spatial distribution of EBs is not only in- 
homogeneous; in several cases one may calcxilate 
reliably a fractal dimension of the EB clusters. 
The fractal behavior of EB clustering occurs for 
any given background definition, provided that the 
statistics are sufficient. Although different defini- 
tions of the background give rise to different esti- 
mations of the fractal dimension, the value of the 
fractal dimension does not change significantly in 
the course of time for a given background. For a 
background contrast Ic 0 = 0.08, the mean fractal 
dimension of EB clusters is ~ 1.4. The relative 
constancy of the fractal dimension suggests that 
the clustering preferences of EBs do not change 
during our observing interval. This is because the 
prolific areas in EB production (the supergranular 
boundaries and the vicinity of the leading spot), 
or the flux emergence process in general, have life- 
times longer than our observing time ( /■v 3.5 h). 

We defined and calculated typical parameters 
of EBs, such as EB areas, maximum and total in- 
tegrated contrasts. To account for unresolved ef- 
fects and the finite spatial resolution, we applied 
a gaussian correction. Independently of the appli- 
cation of this correction, the differential distribu- 
tion functions of all typical EB parameters exhibit 
a robust power-law shape. Moreover, all EB pa- 
rameters are non-linearly interrelated, with corre- 
lation diagrams obeying power laws. This result 
is shown for the first time in an EB study, but it is 
well documented in studies of explosive events in 
the EUV and soft X-rays, as well as in studies of 
hard X-ray flares (Crosby, Aschwanden, & Dennis 
1993; Bromund, McTieman, & Kane 1995; Crosby 
1996; Georgoulis, Vilmer, & Crosby 2001). 

In terms of energetics, EBs have typical sub- 
flaring energies. To account for the calculated ra- 
diative losses and kinematic effects, the peak rate 
of energy release in EBs should vary in the range 
(4 x 10 24 ,10 26 ) erg s -1 , with the total released 
energy in the range (1.8 x 10 27 , 3 x 10 28 ) erg , as- 
suming a typical EB duration of ~ 10 min . We 
have shown that EBs appear in at least two "wave- 
lengths, namely, the wings of the Ha line and the 
1600 A continuum, indicating a vertical extent of 
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at least 100 km. If the bulk of the released en- 
ergy is thermal, then the EB volume is heated at 
a temperature ~ 2 x 10 3 K higher than the ambi- 
ent temperature in the low chromosphere and re- 
sponds to an instantaneous heat deposition with 
a cooling timescale of the order of a few s. The 
distribution function of EB energies obeys a power 
law with index ~ —2.1. While power laws in the 
distribution function of EB energies should be ex- 
pected in any case, the power-law index is sensi- 
tive to the assumptions under which energies axe 
estimated (constant EB lifetime and the gaussian 
correction). 

For decades, EBs have been a puzzling aspect of 
the activity in the low chromosphere. Two ques- 
tions concerning EBs, namely (i) their triggering 
mechanism and (ii) their role on chromospheric 
dynamics and heating, are still elusive. We have 
shown (§ 7) that the locations of EB triggering 
are, quite likely, parts of QSLs in the low chro- 
mosphere. Turbulence plays a key role in the 
evolution of QSLs and the free magnetic energy 
cascades down to small scales, to be released by 
stochastic magnetic reconnection. This intrinsic 
self-similarity fits nicely with the power laws un- 
covered in every aspect of EB statistics and with 
the fractal properties of EB clustering. The typ- 
ical peak EB energies axe similar to the energies 
of the theoretical nanoflares, postulated by Parker 
(1988). Although dependent on various assump- 
tions, the power-law index of EB energies (~ —2.1) 
suggests that EBs contribute significantly to the 
heating of the low chromosphere in the active re- 
gion. If the power-law index is smaller than -2, 
then it can be readily shown that the bulk of 
the integrated energy content of the distribution 
functions stems from the low-energy part (Hudson 
1991). Large numbers of un-detected EBs, there- 
fore, may release substantial amounts of energy. 
It is very likely that our distribution functions fail 
to include the weaker EBs, since typical EB sizes 
axe close to our spatial resolution and numerous 
EBs are hardly detectable. Since EBs occur close 
to the footpoints of flux ropes, it is likely that 
they heat significantly the low chromosphere and 
trigger disturbances that propagate along the field 
lines, thus contributing to a coronal heating pro- 
cess which originates from low-altitude magnetic 
reconnection. 


The idea that EBs are flare-like events is not 
new; it was first proposed in early studies by Sev- 
erny and co-workers (Severny 1968 and references 
therein). Our work, based on a large, seeing-free 
EB sample and high-resolution photospheric ob- 
servations, suggests that magnetic reconnection 
may indeed be a viable triggering mechanism of 
EBs. Several statistical properties of EBs axe sim- 
ilar to statistical properties of EUV and soft X-ray 
BPs in the transition region and the low corona 
(Shimizu 1995; Krucker & Benz 1998; Berghmans, 
Clette, & Moses 1998; Parnell & Jupp 2000; As- 
chwanden et al. 2000; Aletti et al. 2000). Several 
of our conclusions, therefore, could be proved valid 
for these transients, higher in the atmosphere. 

EBs are familiar to experienced solar observers 
but their significance and impact in various as- 
pects of chromospheric dynamics has probably 
been overlooked. A significant effort needs to be 
focused on the study of EBs. Attempts to in- 
terpret EBs by means of phenomenological, or 
MHD/radiative transfer, models, axe highly de- 
sired. If it is shown that EBs are signatures of 
low-altitude reconnection and a viable of 

chromospheric heating, this would be an impor- 
tant advance. Such a result would explain why 
emerging flux regions, typically showing swarms 
of EBs, appear brighter than fully developed _ARs, 
with far fewer EBs and a decaying magnetic evo- 
lution. 
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A. Gaussian correction for finite spatial resolution and unresolved effects in EBs 



Fig. 13 — Application of the gaussian correction for an EB at a given time (frame): The EB contrast 
famts accordmgiy away from the pixel with the maximum contrast at that time. For a background contrast 

,V/ T 7° measure 311 EB area So, which results in a length l 0 , calculated under the assumption of an 
elliptical EB shape with the length 40% larger than the width. The corrected EB length is l m , from which 
we calculate the corrected EB area S m . 


Let us assume that the contrast 1(1) of an EB at any given frame faints as a gaussian function of the 
distance l from the pixel of maximum contrast J m , that is, 


!( l ) = Imexp[- 1 —) 


(Al) 


with b being the width of the gaussian (Fig. 13). 

For a background contrast J Co we measure an EB surface S 0 , that is, an EB length l 0 ~ 1.3(S 0 /-x) 1/2 
assuming that the EB exhibits an elliptic shape with a length l 0 40% larger than the width w 0 . This 
assumption is supported by our results (Fig. 3). Equation (Al), therefore, becomes 


I(S) ~ J m ea:p[— 0.54^-] 


(A2) 


The corrected length of the EB is l m ~ 1.3(5 m /7r) 1 / 2 , with S m being the corrected EB area. Setting the 
background contrast equal to unity and assuming (Fig. 13) that I(S) = 1 for / = l m /2 ~ 0.37Sm /2 , we obtain 


b 2 ~ 


0.375„ 

lnl m 


Substituting equation (A3) into equation (A2) we obtain 


(A3) 


7(5) ~ 4-i-46(s/s m ) 


(A4) 


Moreover, for I(S) — I Co we obtain l = l 0 /2. Substituting this information in equation (A2), one may 
calculate the corrected area S m of the EB at a given frame as a function of the measured area So of the 
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EB at that frame, the peak EB contrast I m at that frame, and the prescribed background contrast I c as 
follows: 


Srn — Sq 


logI„ 


log(l m /lc 0 ) 


(A5) 


The corrected, spatially integrated contrast I s (t), at a given time t (a given frame), is provided by the 
integral 

Is{t) = 2 J 0 1{l)dl (A6) 

Substituting equation (Al) into equation (A6), we obtain 


Is(t) = 2 I m b f 

Jo 


/m/26 


e “ du ; u = \ 

b 


which transforms to 

or, in terms of the corrected area S m , 




I s (t)~bI m VZerf(?l^L.) 

Zb 


(A7) 

(AS) 

(A9) 


with erf (x) = ^ f* e u du being the error function of the quantity x. The coefficient b is given by equation 
(A3). 
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